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ABSTRACT 


(Distribution Limitation Statement No. 3) 


In the process of calculating opacities which are needed as input data for 
radiation transport problems, the absorption cross sections for various species 
has to be known. One factor in the cross section is the vibrational transition 
probability, called the Franck-Condon factor (FCF). This report presents input 
data needed to calculate these from Rydberg-Klein potentials as well as tables 
of Franck-Condon factors for species of atmospheric interest. 
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SECTION I 


INTRODUCTION 


For any given fixed temperature, the intensity distribution of the spectrum 
of a band system in a typical atmospheric diatomic molecule is determined 
primarily by the vibrational transition probability. The intensity distribution 
within a vibrational band is determined by the rotational transition probability. 
This report will be concerned with only the calculation of the vibrational 
transition probabilities called the Franck-Condon factors. These are defined 
by the following overlap integral: , 


2 
Vytyt = ( vy ve cr) (1) 


where yee Vie are the vibrational wave functions to be described in detail 
later and v', v" refer to the upper and lower vibrational quantum nwber, 
respectively. The Franck-Condon factors satisfy the following relations: 


all 2, dyryn tl 
and 


all > qyty tl (2) 


y" 


Since the overlap integral involves the vibrational wavefunctions of tie 
initial and final states of the molecule, the Schrodinger equation must be 


solved for these wavefunctions. 
Hy =» E. y (3) 


where 
%5 is the vibrational wavefunction of level v 


E_ is the energy of vibrational level v 


v 
H is the Hamiltonian operator for the problem 
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The Hamiltonian operator contains a term for the potential energy in which the 
two nuclei move. The correctness of this potential will determine how good the 
vibrational wavefunctions will be. This in turn determines the accuracy of the 
Franck-Condon factors. Hence the better the potential the better the Franck- 


Condon factors. 


There have been a great many empirical potential energy functions developed 
over the past few decades. A study of these and the relative merits of each 
has been done by Steele et al. (Ref. 1). 


Among the most popular potentials in the past for calculating Franck-Condon 


factors has been the Morse potential. 
-ax]? 
V(r) = De E -e ] (4) 


whe re 


bs bn? c2we2y 's 
. 2De 


De is the dissociation energy 
x = (r - re) 


re is the internuclear separation at the minimum of the 
potential well » 


r is the intermuclear separation 


Using this potential one can develop a closed form representation for the wave- 
functions, hence the evaluation of the overlap integrals can be done without 
the need of a large computer. However, this potential is accurate only near 
the bottom of the potential well. Therefore the Franck-Condon factors for the 


upper vibration levels may be inaccurate. 


To calculate better wavefunctions for the upper vibrational levels, one 
needs a better potential. One such method for calculation a better potential 
has been developed by Rydberg (Ref. 2) and Klein (Ref. 3), which makes use of 4 
experimentally measured data. Subsequently, it has been modified in various 
ways by other investigators, notably by Rees (Ref. 4). The general method has . 


become known as the Rydberg-Klein-Rees or R-K-R method. In many instances this 


is a misnomer because the Rees modification is not even used. ¥ 
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The method is based on a WKB approximation and will be outlined below. This 
treatment of the basic equations is based on the derivation in a paper by Van- 
derslice et al. (Ref. 5). 


Rydberg used two expressions from classical mechanics and modified them 
appropriately to obtain the corresponding quantum mechanical expressions. The 
classical equations for the action integral is 


fp, are (2u) g [i = eee" dr (5) 


momentum in r direction 
internuclear separation 
reduced mass of the molecule 
effective potential 


total vibrational-rotational energy given by 


~ (r,) 22 ~ 
Us ——— + —— + V(r) (6) 
2u Qur2 


angular momentum of the rotating molecule 


potential energy of the rotationless molecule 


one then has 
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where 


~ 


K ~ 
Vere = re V(r) 


r 


from which follows equation (7) after solving the P.. and substituting this 
expression in the action integral. | 


The action integral equation (17) is quantized as if it were a harmonic 
oscillator which results in the equation 


h(vt};) = (2u)? § [i - Tore] dr (8) 


This is the first basic equation of Rydberg's method. 


The second basic equation of Rydberg involves the rotational energy of the 
system. Again, classically this is given by 


K 
~ 1 K 1 1 
Eek(i “= S(E peed 2 (9) 
. e ty 7 Ar? V r2P | 
v r 
where 
K is the same as for equation (7) 
Ty is the period of vibration 


Quantization of K the rotational angular momentum leads to 


2 


K = J(J+1) (10) 


8r2y 


with which the rotational energy integral becomes, upon solving equation (15) 
for Pe and substituting this in equation (17), 


h2 1 sf ~ fe dr ih? (+) 
—— BG = ¥ oe eee beet ee (11) 
8r2r ¥2y eff r? 8n2y \r2 . ij 


This is Rydberg's second basic equation. Rydberg solved these graphically. 


The second of these becomes infinite at the classical turning points since 


U = ¥ there. This increases the difficulty of solving the problem graphically 


since large errors can be made in the regions of the classical turning points. 
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Klein defined a new function S(U,K) to aid in finding the classical turning 
points, given by 


I' 
s(U,K) naar J [U - E(1,k))? a1 (12) 


where 
I = h( ves) 


U is the total energy 


22 
ko 


E(1,K) is the vibrational-rotational energy 
I' is the value of I such that U = E(I,k) 


Using this function Klein found that expressions for the classical turning 
points were given by 


3s 
—sfs r 13 
aU ' (ae nin) (13) 
and 
r) 
Sa) -= (14) 
oK min max 
which can be inverted to give 
r, a(f£4 ¢2)- ¢ qs) 
min g 
r. oft e2)4 ¢ (16) 
max g 


The integration and differentiation of Klein's equations must be carried 
out numerically in general. Before the large computers of recent years this 
presented a problem. Rees therefore made the next modification to the Rydberg- 
Klein method, by observing that the integral for $(U,K) could be integrated 
analytically if in the function E(1,K), I is expressable as a quadratic (or 
cubic) form in (vt). Consequently, f,g can be written in closed form. 


Vanderslice et al. capitalized on this observation and have fitted the 
energy E(1,K) with a series of quadratics over the whole energy range of the 
potential function, They then obtained a closed form representation for g,f 


5 
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and consequently, Tain 


LJ -1 
»» ro (u» : 2) 


i=] 


+ (wx); ? (22, - a, (wx), «,) ln | 


a 45 Nya s 
Wy (wx), U, Wy - 2 (wx): Up 


Wy - 2(ux), 4 ud 


ae 
Wi 4(wx), Uy 


Wy is the coefficient of (v+) 
wx, is the coefficient of (v+t';)2 
the coefficient of (v+s) J(J+1) (coupling constant) 


the coefficient of J(J--1) 


Vanderslice et al. also give a simple interpretation of the s(U,K) function 
of Klein, s(U,k) is one half the area between the constant total energy curve, 
0. and the potential curve Voge) (see figure 1). 


Figure 1. Interpretation of the Klein S Function 
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Following Vanderslice, 
dI = tdE (20) 
with 


1. (x)° (e . oa dr (21) 


putting tdE in for dI in Klein's equation for S, 


I 
~ 5 1 ~ ~ I 
S(U,K) ane [U - £(1,k)] at 


one finds, 
f 5 
1 ~ ~ ye ~ 
s “5, (U - £) [( . Yees) | dE 
Ce] 
or 
U 
$= 5dr \ G@-H(E-V,,)i (22) 
v 
eff 


The lower limit of the second integral is Vere since it is evaluated at a 


constant value of r. 
The last integral becomes, after some manipulation, 


r 


max U 
sei dr at - 1, (23) 
2 2 
Thin Vere 


1. METHODS OF VARIOUS INVESTIGATORS 


Many investigators have reported calculating RKR potentials for molecules 
of importance in heated air. Table 1 contains the name of the investigator (or 


investigators) and the general method used. 


] 


1 
{ 
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Table 1 


INVESTIGATORS USING A WKB-TYPE APPROXIMATION 


Flinn, Spindler, Fifer RKR 
Gilmore RK, KDJ 5 
Ginter, Sattino RKR 
Halmann, Laulicht RK 
Jarmain KDJ 
Jain, Shani, Singh RKR 
Krupenie, Weissman RK 
Richards, Barrow RK-RKR 
(hyb rid) 
*Vanderslice et al. RKR 
Zare, Albritton RK 
AFWL RK 


The RK method uses Klein's equations and numerically integrates them in 
some fashion. The RKR method uses the quadratic fit of Rees in some fashion. . 
The KDJ method makes use of a series developed by Jarmain (Ref. 6) based on a 
series expansion of a WKB approximation by Dunham (Ref. 7). With some simpli- 
fying approximations this expansion is the same as a series expansion of Klein's 
equations, which are also based on a WKB approximation. Hurley (Ref. 8) has 
shown the two methods are equivalent and can be deduced from the same basic 


equation, 
2. LIMITATIONS OF RK OR RKR NETHODS 


Ginter and Battino (Ref. 9) have discussed the effects of experimental 
errors and errors due to approximations made in the formulation of equations 
used in the RKR methods on the calculated potential curves. A summary of their 


conclusions will be presented, 


*Co-author with several people on different papers. 
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The accuracy with wnich one can represent the vibrational-rotational energy 
levels of an electronic state of a diatomic molecule by a polynomial E(v,J) 
primarily determines the accuracy of the calculated potential. The errors in 
determining E(v,J) fall into three categories: 


(1) Experimental error in spectroscopic data 


(2) Error caused by fitting the best polynomial E(v,J) to the spectroscopic 
data 


(3) Error caused by fitting the spectroscopic data to a particular form of 
E(v,J) 


ff and AG 


and errors due to using the effective constants directly. From (2) often times 


From (1) there can be errors because of uncertainty in determining BY 


the experimental data will fit with any one of several expansions and do it 
equally well. From (3) if a series of expansions of quadratic form for E(v,J) 
is assumed, then errors may be introduced because of the fitting of local fluc- 
tuations in the experimental data. 


3. ZARE PROGRAM AT AFWL 


Several years ago R. N. Zare made available to the Air Force Weapons Labore- 
tory (AFWL) a set of computer programs needed to construct R-K molecular poten- 
tials for diatomic molecules and to calculate Franck-Condon factors, r-centroids, 
and relative intensities from these potentials. A summary of the methods of 
the programs will be given. 


The molecular potential generation program is one which was written by 
J. Kasper (Ref. 10). The equations, 


I* 
h te 
Sar j U - E(I,K) dI (24) 
2n(2y) 4 ] 
I° 
h ~}5 
g = B. [U - E(I,K)] 7 dI (25) 
2n(2u)? \ vl 


are integrated directly without assuming an analytic form for the integral 
near the upper limit. 
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For the energy E(I,K) of the level of interest either the conventional 


expansion 


E(v,J) = witvek) - wox (vt)? + wy, (vt) 3 


1.5% 5 
+ Woz, (vet's) + wet, (vrs) Ye eo 


+B J(J+1) - D, J2(541)27 +... (26) 


or the CB obtained directly from the spectroscopic data can be used. The 
variables I,K are of course related to V,J by 


I = h(v+}) (27) 


K = J(J+1) (28) 


2n27y 


From the f,g the values of Tain and Trax ate calculated for each vibrational 
level. These are then used as input data to the Franck-Condon factor program. 
To calculate the wavefunctions needed in the Franck-Condon factor program, the 
radial Schrodinger equation is solved (for the wavefunctions) by a numerical 
method developed by Cooley (Ref. 11). The one-dimensional second order 


Schrodinger equation, 


d2 
ore [ew = Vege(Fo3)| y=0 (29) 


where 


E(v,J) is the energy of the particular level, 


h2 


Voge (tod) = V(r) + ; 
Tur 


J(J+1) (30) 


is the energy of rotationless level corrected for rotation, and y, the 


vibration-rotation wavefunction, is changed to a finite difference equation. 
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The finite difference equation is solved iteratively for E(v,J) the eigen- 
value and y the eigenfunction, until the E(v,J) and the previous E(v,J) differ 
by an arbitrarily small amount. The wavefunctions thus obtained are used in 
the evaluation of the square of the overlap integral which is the Franck--Condon 


factor. 
dyty * Sv yey Yer gne dr? (31) 
They are also used in evaluating the r-centroids. 


y rv, dr 
banae LAS dad | (32) 
Svyr vyn de 


¥. (33) 


The integrals are all evaluated using a Simpson's rule numerical integration 


procedure. The relative intensities scaled to ten are also calculated. 


The programs provided by Dr. Zare have been adapted for use on the AFWL CDC 
6600 computer. Since this computer carries 15 decimal numbers per word, the 
double precision subroutines were removed. The (v',v") array size for the 
computed Franck-Condon factor was increased to (20,23). Some of the input 
statements were modified and output statements were added. Aside from these 
changes the program used at AFWL is the same as the one developed and used by 
Dr. Zare. 
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SECTION II 


MOLECULAR BAND SYSTEMS CONSIDERED 


The motivation for calculating Franck-Cordon factors at AFWL is their use 
in calculating transition strengths of diatomic molecules. These strengths 
are used in part to calculate absorption coefficients which in turn are fre- 
quency averaged over prechosen frequency groups to give opacities for air. 
These opacities are desired over a large temperature density range. When deter- 
mining which molecules are candidates for inclusion in the opacity calculation, 
one must take into account the chemistry that goes on in this heated air. 
After picking the potential molecules one looks at the various transitions 
allowed in the molecules and the relative strengths and positions (in energy) 
of the transitions. Then from these one arrives at a list of important elec- | 
tronic band systems of diatomic molecules in heated air. 


The following table contains information on all transitions considered in 
the present work. The transition involved, common name of the band system, 


the geource of the potential, and array size calculated are listed. : 


Table 2 


e 
es Sa Sane 


TRANSITIONS INCLUDED 


Source of 
potential 
State Array size | 
States involved Name Upper Lower (viv") 
CN B2n*-cn x?5* CN violet 13 13 (18,18) 
CN AP ny y-€N X757 CN red 13 13 (18,18) 
CN AZ ny s2-CN X22" CN red 13 13 (18,18) 
Co Aln-co x!5* CO 4th positive 16 16 (20,24) 
Np B?n oN, Aste Ist positive 18 18 (17,13) 
N> C3n -Np B35, 2nd positive 18 18 (4,17) - 4 
| 
No bt lr oN, KID Birge-Hopfield #1 33 19 (20,21) ; 
b' lr -N2 xIye Birge-Hopfield #2 33 19 (24,21) 
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Table 2 (cont'd) 


Source of 
potential 
State Array size 
States involved Name Upper Lower (v'v") 
. Nz A2n -N> oe Meinel 330233 (20,20) 
nt B2r"-N3 oe Ist negative 33 33 (20,20) 
Nz C257-N3 x2r° 2nd negative 33033 (10 ,20) 
No A2z*-NO X2m) 2 y (gamma) 22.23 (5,23) 
NO A2r*-NO X2n3 72 y (gamma) 22 «28 (5,23) 
NO B2n-NO X2nj 72 B (beta) 22. «23 (19 ,23) 
NO B2n-NO X*13 2 8 (beta) 22 23 (19,23) 
NO C?n-NO X?m) 2 6 (delta) 22 23 (4,23) 
NO C?n-NO X?m3 2 § (delta) 22 23 (4,23) 
NO D2z*-NO x2, /2 ¢ (epsilon) 33.23 (7,23) 
No D?E*-No X2n3 /2 ¢ (epeiion) 33023 (7,23) 
I NO B'7A,-NO X?m 7, 8° (beta prime) 22 23 (5,23) 
. NO B'24,-NO X?n5 2 8° (beta prime) 22 23 (5,23) 
NO E25*-No X20 /2 y' (gamma prime) 22 23 (5,23) 
NO E2E"=NO X2n3 79 y' (gamma prine) 22 (5,23) 
NO G2E-NO X?mq 72 Lagergmist-Miescher 33 23 (7,23) 
NO C2n-No A2z* Heath 22.22 (4,5) 
NO D2Z*-No A2r* Feast 1 22 22 (7,5) 
NO E2r*-No a2z* Feagt 2 22.22 (5,5) 
NO E2z*-NO C2n 22 22 (5,4) 
NO E257-No D25* Feast-Heath 22 33 (5,7) 
0, B35, -02 x3Ee Sch umann-Runge 33 33 (20,23) 


+ 
No-No 


KIT Np xr Ist 33 33 (20,20) 


AFWL-TR- 70-108 


Table 2 (cont'd) 


Source of 
potential 
State 
States involved Name Upper Lower 
N> KID} A?n 2nd 3333 
No KID N? B?rr 3rd 33033 
N> KID Ny er 4th 3333 
0,-05 

0> x315-07 X?r Ist 33 33 
0, x35 -07 a‘n 2nd 33.33 
0, X35 -07 A2y 3rd 33033 
0» x35 -07 ber, 4th 33 33 
0, x30 -07 Ee 5th 3333 
0, x325-07 chre 6th 3333 


1. CN (CYANOGEN) BAND SYSTEMS 


Array size 


(v'v") 


(20 ,20) 
(20,20) 


(10 ,20) 


(14,23) 
(10 ,23) 
(10,23) 
(10 ,23) 
(10,23) 
(10 ,23) 


The potentials used to calculate Franck-Condon factors were taken from the 


paper by Fallon, Vanderslice, and Cloney (Ref. 13). Tables 3 through 6 contain 


the data for the potential energy curves for the various CN states used. 


Table 7 contains other molecular data needed by the FRANKON program. 
Table 3% 


POTENTIAL ENERGY FOR THE CN x27 STATE 


v Wem!) Fin Fmex (ev) 
0 1031.1 1.1237 1.2244 0.1278 
1 3073.6 1.0923 1.1676 0.3811 
2 5089.7 1.0722 1.2997 0.6310 
3 7079.5 1.9566 1.3274 0.8777 
4 9042.8 1.0438 1.3525 1.1211 
5 10980.0 1.0327 1.3760 1.3613 


6 12890.0 1.0230 1.3984 1.5981 


AFWL-TR- 70-108 . | 
Table 3* (cont'd) 
v V(em™!) ain * max Viev) | 
. 7 14775 .0 1,0143 1.4198 1.8317 | 
8 166 32.0 1.0063 1.4406 2.0620 
9 18463.0 0.9991 1.4610 2.2890 
10 20267.0 0.9924 1. 4809 2.5127 
11 22045.0 0.9861 1.5005 2.7330 
12 23797.0 0.9803 1.5200 2.0503 
13 25520.0 0.9749 1.5393 3.1639 ) 
14 27217.0 0.9696 1.5583 3.3743 
15 28887.0 0.9648 1.5774 3.5813 
16 30530.0 0.9602 1.5964 3.7850 
17 32146.0 0.9558 1.6153 3.9853 
18 33734.0 0.9516 1.6343 4.1823 
*Reference 13. 
Table 4* 


POTENTIAL ENERGY FOR THE CN A?13/2 STATE 


y v(en"!) — Fmtn Tmax Ye) 

0 903.99 1.1827 1.2900 0.1121 

! 1 2691.1 1.1482 1.3354 0.3336 
| 2 4452.9 1.1269 1.3701 0.5521 
3 6189.4 1.1105 1.4001 0.7673 

4 7900.4 1.0969 1.4273 0.9795 

; 5 9585.0 1.0852 1.4529 1.1884 

6 11245.0 1.0749 1.4773 1.3942 

‘ 7 12880.0 1.0657 1.5007 1.5968 
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Table 4* (cont'd) 


y Veen!) Frntin ‘max = Ycev) 

8 14487.0 1.0573 1.52% 1. 7960 : 
9 16070.0 1.0477 1.5459 1.9923 

10 17627.0 1.0526 1.5679 2.1853 , 
11 19158.0 1.0360 1.5896 2.3751 

12 20663.0 1.0298 1.6111 2.5617 

13 22142.0 1.0240 1.6325 2.7451 

14 23594.0 1.0186 1.6538 2.9254 

15 25021.0 1.0134 1.6750 3.1020 

16 26422.0 1.0083 1.6960 3.2757 

17 27795.0 1.0038 1.7175 3.4460 

18 29142.0 0.9994 1, 7388 3.6130 
*Reference 13, 

Table 5* 


POTENTIAL ENERGY FOR THE CN A? my 72 STATE 


y Vom!) Fin Fmax Yew) 
0 903.99 1.1816 1.2891 0.1121 
1 2691.1 1.1484 1.3356 0.3336 
2 4452.9 1.1272 1.3704 0.5521 
3 6189.4 1.1109 1.4006 0.7673 
4 7900.4 1 5975 1.4280 0.9795 
5 9585.9 1.0860 1.4527 1.1884 
6 11245.0 1.0759 1.4783 1.3942 


i 12880.0 1.0668 1.5019 1.5968 
8 14487.0 1.0587 1.5249 1. 7960 i 


9 16070.0 1.0512 1.5475 1.9923 
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Table 5* (cont'd) 
y V¢em"?) Fun Fmax — “(ev) 
10 17627.0 1.0443 1.5696 2.1853 
ll 19158.0 1.0379 1.4916 2.3751 
12 2066 3.0 1.0319 1.6133 2.5617 
13 22142.0 1.0263 1.6348 2.7451 
14 23594.0 1.0211 1.6563 2.9254 
15 25021.0 1.0161 1.6778 3.1020 
16 26422 .0 1.0115 1.6992 . 3.2757 1 
17 27795.0 1.0070 1.7207 3.4460 
18 29142.0 1.0028 1.7423 3.6130 
*Reference 13. 
Table 6* 


POTENTIAL ENERGY FOR THE CN B2I* STATE 


y Vien!) Fin Tmax (ev) 
0 1077.0 1.104 1.203 0.1335 
1 3200.0 1.074 1,246 0.3968 
2 5284.0 1.055 1.279 0.6551 
3 7325.0 1.040 1.307 0.9081 
4 9321.0 1.027 1.333 1.1556 
5 11269 .0 1.016 1.358 1.3971 
6 13167.0 1,007 1.383 1.6324 
7 15011.0 0.999 1.407 1,811 

3 8 16799.0 0.992 1.432 2.0827 
9 18527.0 0.986 1.457 2.2969 

10 20192.0 0.980 1.483 2.5034 
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Table 6* (cont'd) 


y Vem"! “ain 
11 21795.0 0.974 
12 23338.0 0.970 
13 24821.0 0.965 
14 26249 .0 0.961 
15 27623.0 0.957 
16 28948.0 0.953 
17 30228.0 0.950 
18 31464.0 0.946 
*Reference 13. 

Table 7 


ev) 

2.7021 
2.8933 
3.0773 
3.2542 
3.4247 
3.5889 
3.7475 
3.9008 


ADDITIONAL DATA NEEDED FOR THE CN MOLECULE 


State De be 

x25 61525.88 1.8996 
A?n3/2 52279.86 = 1.7165 
A?r, 72 $2279.86 = 1.7165 


Ber 54774. 34 1.9701 


2. CO (CARBON MONOXIDE) BAND SYSTEMS 


ae 
0.01735 
0.01746 
0.01746 


0.02215 


¥ 
0 


° 


Te 
0.0 
9222.04 
9272.62 
25752.4 - 


The only band of CO considered to be of possible importance in absorption 


in heated air is the CO 4th positive system. 


The potentials used in calculating 


the Franck-Condon factors were taken from the paper by Krupenie and Weisman 


(Ref. 16). Tables 8 and 9 contain the data for the potential curves of the 


lower and upper states. Table 10 contains other molecular data needed by the 


program. 
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Table 8* 


POTENTIAL ENERGY FOR THE Co x!z* sTaTE 


y V(cm™!) mtn ‘max = (ev) 
0 1081.590 1.083 1.179 0.1341 | 
1 3224. 864 1.053 1.220 0.3998 | 
2 5341.654 1.034 1.250 0.6623 
3 7432 .026 1.020 1.276 0.9214 
4 9496 .948 1.008 1.300 1.1773 
5 11533. 76 0.997 1,322 . 1.4300 
6 13545.20 0.988 1.343 1.6794 
7 15530.64 0.980 1.363 1.9255 
8 17489 .9 0.972 1.383 2.1684 
9 19423.5 0.965 1.402 2.4081 
10 21330.7 0.959 1.420 2.6446 
ll 23212.4 0.953 1.438 . 2.8779 
12 25068.5 0.947 1.456 3.1080 
13 26898.4 0.942 1.474 3.3349 
14 28703.2 0.937 1.492 3.5587 
15 30482.7 0.932 1.509 3.7793 
16 32235.4 0.928 1.526 3.9966 
17 33963.9 0.923 1,544 4.2109 
18 35666.8 0.919 1.561 4.4220 
19 37344 .4 0.916 1.579 4.6300 
20 38997.5 0.912 1.596 4.9349 
21 40625.7 0.908 1.614 5.0368 
22 42208.0 0.905 1.631 5.230 
23 43776 .0 0.901 1.649 5.4274 


*Reference 16. 
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Table 9% 


POTENTIAL ENERGY FOR THE CO A*n STATE 


y V¢cm"!)  Fmin Fux = V(ev) 
0 753.49 1.183 1,297 0.0934 

1 2202.3 «1.148 = 14348 (0.2780 
2 3685.1 1.126 1.387 0.4569 
3 5097.9 1.110 1.422 0.6320 
4 6476.1 1,096 1.454 0.8029 
5 7818.2 1.085 1.484 + —«-0.9693 
6 9125.0 1.075 1.516 1.1313 
7 1041.8 1,066 1.542 1.2896 
8 11641.0 1.058 = 1,571 1.4433 
9 12846.9 1.051 1.599 ~—s:1.5928 
| 10 14018.8 1.046 1.627. 1.7382 
ul 15155.6 1.038 1.656 ~—S>_ 1.8790 
12 16262.0 1.031 1.685 2.0162 
13 17333.0 «1.02772, 1490 
14 18371.0 1,022 1.743 -2.2777 
15 1936.0 1.017 1.7732, 4010 
16 20330.0 1.013 1.804 2.5205 
17 21257.0 1,009 1.837 2.6355 
18 22.38.0 1.005 1.871 2.7447 
19 22970.0 1.001 1.090 —«-2, 8478 
20 23762.0 0.996 1.051 2.9460 


*Reference 16. 
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Table 10* 


2] ———<—— TD 


ADDITIONAL DATA NEEDED FOR THE CO MOLECULE 


(u = 6.85841) 


State De Be ae ye ‘Te 
xist = .g9597.4 1.93127 0.617513 2.960E-6 0.0 
Aln 24522.8 1.226 0.018 0.0 133740.0 


*Gilmore, DASA Handbook 1917 old copy Volume I tables. 


3. No (NITROGEN) BAND SYSTEMS 

The potentials used to calculate the Franck-Condon factors for the nitrogen 
systems came from several sources. The nr, Aart, Bens and C3n states 
contained in tables 11 through 14 were taken from the paper by Benesch et al. 
(Ref. 18). The potentials used for the bin and a states were generated 
at AFWL with program TURNGPT and are listed in tables 15 and 16, respectively. 
The data for the bin state are from Lofthus (Ref. 19) and the data for the 


birt state are from Wilkinson (Ref. 20). These data are listed in tables 17 


through 18, 
Table 12% 
POTENTIAL ENERGY FOR THE Np xine STATE 
v V(cm7!) Tain (A) * max (A) View) 
0 1175.5 1.0559 1.145, 0.14573 
1 3505.2 1.026, 1.184, 0.42456 
2 5806.5 1.008; 1.213 0.71987 
3 8079.2 0.9945 1.2377 1.00163 
| 4 10323.3 0.982, 1.2605 1.27985 
aie 5 12538.8 0.9725, 1.280, 1.55452 
6 14725.4 0.963, 1.3007 1.82561 
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Table 11* (cont'd) 


v veem!) Fain’) Fax) Vey) 

7 16883.1 0.955, 1.3197 ~—«-2.09311 

8 19011.8 0.948, 1.338) 2.35702 

9 21111.5 0.9417 1.356, 2.61733 
10 23182.0 0.935; 1.373, «2.87403 

uu 25223.3 0.929, 1.391; +3270 

12 27235.3 0.924, 1.408) 3.37654 

13 29218.0 0.919, 1.425, 3.62235 
1 31171.2 0.914, 1.4429: 86.450 
15 33094.9 0.910, 1.458, 4.10300 
16 34989.0 0.905, 1.4753 4. 33782 
17 853.5 0.9017 1.492; 4 56.898 
18 38688.3 0.897 1.508, 4. 79645 
19 40493.4 0.894) 1.5255 «5.02024 
20 42268.6 0.890, 1.5423 5.24032 
21 46014.1 0.8879 «1.589, 5.45672 


*Reference 18. 


Table 12* 


POTENTIAL ENERGY FOR THE N, Aare STATE 


v V(om7!) Tain (A) Tmax A) Viev) 
0 726.80 1.2329 1.3482 0.0901 
1 2159.67 1.1979 1.3992 0.2678 
2 3564.87 1.1758 1.4378 0.4420 : 
> 4942.36 1.1589 1.4715 0.6127 
4 6292.03 1.1450 1.5025 0.7801 ‘ 
5 7613.75 1.1332 1.5319 0.9439 
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Teble 12% (cont'd) 


v V(cm™!) 
6 8907.33 
7 10172.54 
8 1140910 
9 12616 .69 
10 13794.95 
11 14943.47 
12 16061. 80 
13 1714944 


*Reference 18, 


POTENTIAL ENERGY FOR THE N, 


y V(cm™! 
0 863.16 
1 2568.14 
2 4244.46 
3 5891.92 
4 7510.38 
5 9099.74 
6 10659.92 
7 12190. 86 
8 13692.50 
| 9 15164. 76 
10 16607.54 
fo: 11 18020. 70 
12 19404.04 


Faint) 
1.1229 
1.1138 
1.1055 
1.0980 
1.092, 
1.0845 
1.078, 


1.0732 


Table 13* 


Tain OA) 
1.1631 
1.1306 
1.1099 
1.0940 
1.0809 
1.0697 
1.0599 
1.0511 
1.0432 
1.0360 
1.0294 
1.0233 
1.0176 
23 


Tmax (A) 


1.5602 


1.5877 
1.6149 
1.6417 
1,668. 
1.695, 
1,722, 


B3n STATE 
Tax (A) 


1.2689 
1.3150 
1.3497 
1.3798 
1.4074 
1.4334 
1.4583 
1.4825 
1.5061 
1.5294 
1.5524 
1.5753 


1.5981 


“(ev) 

1.1043 
1.2612 
1.4145 
1.5642 
1.7103 


1.8526 


1.9913 


2.126) 


Yew) 
0.10702 
0.31840 
0.52623 
0.73049 
0.93114 
1.12819 
1.32163 
1.51144 
1,69 761 
1.88014 
2.05902 
2.23422 


2.40573 
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Table 13* (cont'd) 


v V(em7!) Fain (A) Taax A) Viev) 

13 20757.31 1.012, 1.621 2.57351 

14 22080.16 1.007, 1.644 2.73752 

15 23372.16 1.0027 1.6675 2.89770 : 
16 24632.76 0.998, 1.690, 3.05399 

17 25861.32 0.9949 1.7143 3.20631 


*Reference 18. 


Table 14% 


POTENTIAL ENERGY FOR THE No C3n STATE 
e 


° 
¥ V(cm7!) Tain“) = “max (A) View) 
0 1016.71 1.1030 1.2005 0.12605 
1 3011.11 1.0730 1.2438 0.37332 
2 4951.90 1.0540 1.2771 0.61394 
3 6825.93 1.0395 1.3076 0.84629 ; 
4 8607.21 1.0282 1, 3388 1.06713 
*Reference 18. 
Table 15* 

POTENTIAL ENERGY OF THE N> bin STATE 
y vem?) Fin Fmax (ev) 
0 412.5 1.23007 1.38773 s:0..0511 
1 1114.5 1.16992 1.45950 0.1381 
2 1779.5 1.13007. 1.51097 =: 0.2206 
3 2407.5 1.09562 1.55852 ° 0.2984 . 
4 2992.5 1.06450 1.60407 0.3710 
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Table 15* (cont'd) 


y Vem"! Fain = Tmax = “(ev) 
‘ 5 3546.5 1.03665 2.64718 0.4396 

6 4062.5 1.00840 1.69324 0.5036 

7 4542.4 0.98508 1.73570 0.5631 


*From TURNGPT using Lofthus data. 


Table 16% 


POTENTIAL ENERGY OF THE N, b'!r’ STATE 


v V(cm7! Fnin "max View) 
0 374.6 1.36740 1.52792 0.0464 
1116.6 1.31442 1.59394 0.1384 
2 1848.9 1.389 36 1.65220 0.2292 
3 2571.6 1.25162 1.68331 0.3188 
4 3284.7 1.22812 1.720% 0.4072 
5 3988.2 1.20746 1.75474 0.4944 
6 4682.0 1.18887 1.78724 0.5804 


*From TURNGPT using Wilkinson data. 


Table 17 


ADDITIONAL DATA NEEDED FOR THE No MOLECULE 


State De Be s xe Te 
xine 78717.0 1.9987 0.0178 0.0 0.0 
| ; Ajre 28499 .0 1.4545 0.01794 -9.2 x 107° 50218.0 
B3n 38306 ..0 1.6375 0.01794 -7.4 x 107° 59639.0 
| : Cn 28643.0 1.8259 0.0197 0.0 69 302.0 
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Table 17 (cont'd) 
De Be ae Te 


55655 .0* 1.154 0.0048 104473.0 


bilze 55655.0% 102288.0 


*Estimated value. 


Table 18 
MOLECULAR DATA NEEDED FOR TURNGPT (\ prise ) 
] 
State We Wexe Weye Weze Be ae 


prize 751.64 4,82 0 0 1.154 0.0048 


Table 19 


MOLECULAR DATA NEEDED FOR TURNGPT (: os,) 


yoy By 
0 0 1.41 
l 1.41 
2 1.40 


3 1.41 


4, Nd (NITROGEN MOLECULAR ION) BAND SYSTEMS 

The potentials used to generate the Franck-Condon factors for the No band 
systems were calculated with program TURNGPT. The A’n and i states were 
generated using spectroscopic constants given by Gilmore (Ref. 21). These are 
given in table 24, The x and Best states were generated by using spectro- 
scopic data from Lofthus (Ref. 19). These data are listed in tables 25 and 26. 
The BY for the missing levels was obtained by graphical interpolation. The 


26 


AFWL-TR- 70-108 


generated potentials are listed in tables 20 through 23. 


Table 20* 


POTENTIAL ENERGY FOR THE Nt xen STATE 


I< 


o 


20 


Véoml) 
1099.4 
3274.2 
5416.4 
7525.8 
9602.2 

11645.4 

13654.4 

1563.1 

17570.0 

19477.0 

21349 .0 

23185.9 

24987.0 

26751.8 

28479 .9 

30171.0 

31824.5 

33439.9 

35016.7 

36554.0 


38051.7 


Tain 


1.07794 
1,04341 
1.02546 
1.01122 
0.99944 
0.97947 
0.97947 
0.97184 
0.96441 
0.95803 
0.95185 
0.94604 
0.94066 
0.93552 
0.93092 
0.92638 
0.92239 
0.91853 
0.91482 
0.91152 


0.90810 


*From Lofthus data. 


r 
max 


1.1664 

1.20672 
1.23763 
1.26388 
1.28783 
1.20991 
1.33082 
1.35192 
1.37203 
1.39205 
1.41160 
1.43092 
1.45022 
1.46943 
1.48894 
1.50835 
1.52823 
1.54825 
1.56848 
1.58931 
1.61014 
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Table 21* 


POTENTIAL ENERGY FOR THE N> A’n | STATE 


V(em=!) Tain Tmax Veev) 


947.96 1.02650 1.22743 0.1175 
2821.38 1.09570 1.27158 0.3498 
4664.80 1.07615 1.30488 0.5783 
6478.22 1.06118 1.3338 0.8031 
8271.64 1.04887 1.38479 1.0242 

10015.06 1.03838 1.38479 . 1.2416 
11738.48 1.02917 1.40847 1.4553 
13431.90 1.02097 1.43137 1.6652 
15095. 32 1.01357 1.45371 1.8714 
16728.74 1.00681 1.47563 2.0739 
18332.16 1.00061 1.40726 2.2727 
19905.58 0.99487 1.51867 2.4678 
21449 .00 0.98954 1.53995 2.6791 
22962 .42 0.98455 1.56116 2.8468 
24445 . 84 0.97987 1.58235 3.0307 
25897.26 0.97546 1.60356 3.2109 
27322 .68 0.97130 1.62484 3.3880 
28716.10 0.967356 1.64622 2.5601 
30079 .52 0.96361 1.66774 3.7291 
31412.94 0.96003 1.68943 3.8944 


32716. 36 0.95662 1.71132 4.0560 


*From Gilmore's constants. 
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Table 22* 


POTENTIAL ENERGY FOR THE Nt Br" STATE 


r 


V(cm7!) nin Fax Vey) 


1197.15 1.03331 1.12304 0.1484 

3588.66 1.00546 1.16196 0.4424 

5887.46 0.98795 1.19199 0.7299 

9147.86 0.97409 1.21811 1.0101 
10344,26 0.96323 1.24330 1.2824 
12467.06 0.95501 1.26903 1.5456 
14508.06 0.94720 1.27401 1. 7986 
16459. 16 0.94043 1.32000 2.0405 
18291. 36 0.93404 1.34765 2.2677 
20018. 26 0.92856 1.37643 2.4818 
21614,96 0.92315 1.40798 2.6797 
23094. 86 0.91867 1.44118 2.8632 
24466.26 0.91496 1.117632 3.0332 
25742 .56 0.91145 1.51208 3.1914 
26938. 86 0.90846 1.54832 3.3397 
28065. 46 0.90544 1.58452 3.4794 
29132 .56 0.90303 1.62095 3.6117 
30148.06 0.90104 1.65763 3. 7376 
31114.06 0.89893 1.69439 3.8574 
320 36 .06 0. 89687 1.73160 3.9717 


32858.06 0.89012 1.78474 4.0736 


*From Lofthus data. 
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Table 23% 


POTENTIAL ENERGY FOR THE Nt ong? STATE 


y V(em™!)  Fntin ‘max = “(ev) 
0 1034.08 1.21465 1.311246 «0.1282 
1 3085.06 1.18158 1.34967 0.3824 
2 5111.94 1.15962 1.37772 «0.6337 
3 7113.28 1.14238 1.40182 0.8818 
4 9087.64 1.12798 1.42383 1.1266 
5 11033.58 1.11558 1.44462 1.3679 
6 12949.66 1.10471 1.46470 «1.6054 
7 14832.44 1.09508 1.48437 «1, 8391 
8 16686.48 1.08646 1.50386 2.0687 
9 1850434 1.07874 1.52333 2.2941 
10 20286.58 1.07182 1.54294 2.5150 


*From Gilmore's constants 


Table 24 


ADDITIONAL DATA NEEDED FOR THE Nt MOLECULE 


State We Wexe Weye Weze WeTe Te 
re 2207.17 16.146 -2.85 x 10°2 9,2 x 1074 ) 0.0 

Ae 1903. 42 15.0 0.0 0.0 ) 8219.6 

Br" 2418.7 22.53 -6.70 x 10°! 4.0 x 10~2 0 25461.6 

cent 2073 10.97 -2,40 x 107! 0.0 0  64607.5 

State = Be ae ye os re 

xr 1.9258 0.01743 -1.64 x 107" Oo 1.918 

An 1.748 0.02 0.00 0 1,117 

Bry 2.083 0.0183 -1.65 x 1073 0 1.075 : 
cry 15124 0.0017 -7.50 x 107% 0 1.262 
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Table 25 


MOLECULAR DATA NEEDED FOR TURNGPT (c x2 it) 


ca & 


0 1099.4 1.922 
1 3274.2 1.902 
2 5416.4 1.879 
3 7525.8 1. 861 
4 9602.2 1, 841 | 


5 11645.2 1,826 
6 13654.4 1.808 
7 1560.1 1.781 
8 17570.0 1. 766 
9 19477.0 1.740 
10 21349 .1 1.724 
ll 23185.9 1.703 
12 24987.0 1.683 
13 26751.8 1.663 
14 28479 .9 1.639 
15 3171.0 1.620 
16 31824.5 1,593 
17 33439 .9 1.572 
18 35016.7 1.548 
19 36554.0 1.522 


20 38051.7 1.500 
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Table 26 


MOLECULAR DATA NEEDED FOR TURNGPT N> Bert 


Ba Gv By 
0 1197.16 2.073 
1 3568.66 2.049 : 
2 5887.46 2.025 
3 8147.86 2.002 
4 10344.26 1.968 
5 12467.06 1.926 
6 14508.06 1. 896 
7 16459. 16 1. 896 
8 18291. 36 1.81 
9 20018, 26 1.761 
10 21614.06 1.710 
ll 23094, 86 1.653 
12 24466 .26 1.595 1 
13 25742.56 1.545 
14 269 38, 86 1.494 
15 28065. 46 1.452 
16 29132.56 1.404 
17 30148.06 1,361 
18 31114.06 1.323 
19 320 36 .06 1.285 
20 32858.06 1.250 
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5. NO (NITRIC OXIDE) BAND SYSTEMS 


The molecular potentials used to calculate the various NO band system 
Franck-Condon factors came from three sources. The Kem pos X20 3 72 states were 
from Vanderslice (Ref. 23). The art, B2n, C2n, B'2A, E25* states were from 
Vanderslice et al. (Ref. 22). The NO D?* state potential was generated with 
program TURNGPT using molecular constants from Gilmore (Ref. 21). The 
potentials are listed in tables 27 through 34. The spectroscopic constants 
used to generate the D state potential are listed in table 35. 


Table 27* 


POTENTIAL ENERGY FOR THE NO X2m /2 STATE 


V¢em™!) Tmax (A) Tain A) Viev) 
948.5 1.202 1.105 0.1176 
2825.0 1.245 1.074 0.3503 
4673.0 1.276 1.055 0.5795 
6493.0 1,304 1.041 0.8051 
8285.0 1.329 1.028 1,027 
10048.0 1.352 1.018 1.246 
11784.0 1.375 1.009 1,461 
13492.0 1.396 1.000 1.673 
15172.0 1.417 0.993 1.881 
16824.0 1.438 0.986 2.086 
18448.0 1,458 0.974 2.388 
20044.0 1.478 0.980 2.485 
21613.0 1.497 0.974 2.680 
23152.0 1.517 0.968 2.871 
24665.0 1.537 0.963 3.058 
26149.0 1.556 0.958 3.131 
27605.0 1.576 0.954 3.423 


290 34.0 1.596 0.945 3.600 
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Table 27% (cont'd) 
V(cem™!) 
3434.0 
31805.0 
33141.0 
34443.0 
35714.0 


36953.0 


*Reference 23. 


Table 28* 


POTENTIAL ENERGY FOR THE NO X25 


3/2 STATE 


ween?) Fmax(*?) Fain” Veew) 
948.5 1.106 0.1176 
2825.0 1.076 0.3503 
4673.0 1.057 0.5795 
6493.0 1.042 0.8051 
8285.0 1.030 1.027 
10048.0 1.019 1.246 
11784.0 1.010 1,461 
13492.0 1.002 1.673 
15172.0 0.994 1.881 
16824.0 0.987 2.086 
18448.0 0.981 2.288 
20044.0 0.975 2.485 
21613.0 0.969 2.680 


23152.0 0.964 2.871 
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Table 28* (cont'd) 
V(cm7! Tmax (A) 
24663.0 1.538 
26145.0 1.558 
27599 .0 1.577 
29024.0 1.597 
30423.0 1.617 
31790.0 1.637 
33124.0 1.657 
34424.0 1.678 
35695.0 | 1.698 


369 38.0 1.719 


*Reference 23. 


Table 29* 


POTENTIAL ENERGY FOR THE NO AZZ‘ STATE 


view!) Fax) Fain” Yeeyy 
1183.0 1,109 0.1467 
3525.0 1,147 0.4371 
5835.0 1.175 0.7235 
8111.0 1.199 1.006 

10354.0 1.220 1.284 


12564.0 1,242 1.558 


*Reference 23. 
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POTENTIAL ENERGY OF THE NO Bn STATE 
e eo 
r__ (A) Tain 


v V(cm7!) 
0 516.8 
1 1540.0 
2 2549.0 
3 3545.0 
4 4519.0 
5 5492.0 
6 6447.0 
7 7390.0 
8 8323.0 
9 9238.0 
10 10144.0 
1l 11030.0 
12 11911.0 
13 12772.0 
14 13609.0 
15 14463.0 
16 15293.0 
17 16081.0 
18 16859.0 
19 17612.0 


*Reference 23. 


Table 30* 


max 

1.488 1.356 
1.545 1.315 
1.588 1.289 
1.625 1.269 
1.658 1.253 
1.690 1.239 
1.719 1,226 
1.748 1.215 
1.776 1.205 
1.804 1.196 
1.831 1.188 
1.857 1.180 
1,883 1.173 
1.909 1.167 
1.935 1.161 
1.962 1.154 
1.989 1.147 
2.015 1.142 
2.041 1.136 
2.068 1.150 
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Yew 
0.0641 
0.191 
0.316 
0.440 
0.560 


0.681 


0.799 
0.917 
1.032 
1.146 
1,258 
1.368 
1.477 
1.584 
1.688 
1.793 
1.896 
1.994 
2.091 
2.184 
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Table 31* 


POTENTIAL ENERGY OF THE NO C2n STATE 


v V(cm™! Tex) fain (A) Vev) 
0 1194.0 1.109 1.022 0.1481 
1 3556.0 1.148 0.997 0.4409 
2 5922.0 1.178 0.981 0.7343 
3 8192.0 1.204 0.970 1.016 
4 10491.0 1.228 0.961 1.301 
*Reference 23. 
Table 32* 

POTENTIAL ENERGY OF THE NO B'2A STATE 
v V(cm=! Tax (A) Tain A) View) 
0 604.0 1. 368 1.246 0.07493 
1 1789.0 1.424 1.209 0.2219 
2 2943.0 1.466 1.18 0.3649 
3 4062.0 1.503 1.168 0.5037 
4 $177.0 1.539 1.153 0.6420 
5 6212.0 1.570 1.139 0.7703 


*Reference 23. 


37 


AFWL-TR- 70-108 


Table 33* 


POTENTIAL ENERGY OF NO E25* STATE 


v V(em™!) *max cA) Fmin cA) Viev) ; 
0 1183.0 1.112 1.025 0.1467 

1 3525.0 1.150 0.997 0.4371 

2 5835.0 1.177 0.980 0.7235 

3 8114.0 1,202 0.966 1.006 

4 10360.0 1.223 0.955 1.285 

5 12575.0 1,244 0.945 1.559 


*Reference 23. 


Table 34% 


POTENTIAL ENERGY FOR THE NO p2s* STATE 


v V(cm=! ‘min Fmax Vey 

0 1156.336 1.02039 1.10891 0.1433 ; 
1 3438004 0.99320 1.14762 0.4262 
2 5687.031 0.97601 1.17645 0.7050 

3 7918.479 0.97601 1.20099 0.9817 

4 10152.686 0.96318 1.22264 1.2586 

5 12415.274 0.94545 1.24187 1.5392 

6 14737.141 0.93970 1.25881 1.8270 

7 17154.469 0.93579 1.27350 2.1267 


*From TURNGPT using Gilmore's constants 
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Table 35 
MOLECULAR DATA NEEDED FOR TURNGPT (NO D*z) 
We Wexe Weye Weze Be ae ye re 
} | 
2323.9 22.885 0.75 0.22 2.0026 2.175 x 1072 1) 1.062 


6. O 2 (OXYGEN) BAND SYSTEM 


Franck-Condon factors were calculated for only the 0, Schumann-Runge band 
system. The potentials needed, 0, X°r and 0, B35, were calculated with program 
TURNGPT using data from Wallace (Ref. 24). The G, AG, BY were each plotted and 
the best smooth curve was drawn through these. These smoothed values were 
then used as input data to TURNGPT. The B°5 state was found to have a hook 
near the dissociation limit on the inner branch of the potential curve. This 
is similar to the findings of many other investigators, for example, Vanderslice 
et al. (Ref. 25), F. Gilmore (Ref. 26), Richards and Barrow (Ref. 27),* Ginter, 
Battino (Ref. 28). 


A potential which is not single valued causes numerical difficulties when 
one attempts to find solutions for the wave functions needed to calculate 
overlap integrals. To avoid this problem the calculated potential for the B?r 
state was used up to the vibrational level where the hook started to occur. 
First an extrapolation of the form a/r!? wes used. (This procedure is similar 
to that used by Halmann and Laulicht (Ref. 29) who used Vanderslice's potential 
for their calculations.) Using the resulting potential to find the maximum of 
the Shcumann-Runge continuum by the reflection method, one finds that it occurs 
at a much too high energy. An alternate method for constructing a more realistic 
potential in the troublesome region involves using the calculated potential at 
values of vibrational quantum number where one believes the calculation and using 
Evans and Schexander's (Ref. 30) values for the potential in the continuum 
portion of the potential and interpolating through the troublesome region. 
Franck-Condon factors were calculated using the so determined potentials. 


*They claim no anomaly though. 
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Table 36* 


POTENTIAL ENERGY FOR THE 0> x7r) STATE 


V(cm™! 
789.81 
2352.20 
3891.68 
5408.25 
6901.91 
8372.64 
9820.43 
11245.26 
12647.,12 
14025.99 
15381. 88 
16714. 79 
18024.69 
19311.59 
20575.48 
21816. 33 
2334.10 
21228, 74 
25400.22 
26548.54 
27673.62 
28775.15 
29852 .56 


30905.29 


*From TURNGPT 


Tain 
1.15966 
1.12759 
1.10739 
1.09177 
1.07892 
1.06791 
1.05849 
1.04980 
1.04186 
1.03482 
1.02825 
1.02215 
1.01657 
1.01136 
1.00645 
1.00183 
0.99764 
0.99367 
0.98988 
0.98618 
0.98325 
0.97498 
0.97249 


0.97124 


— 
1.26310 
1.30786 
1.34156 
1.37059 
1.39709 
1.42192 
1.44586 
1.46867 
1.49078 
1.51266 
1.53410 
1,55527 
1.56737 
1.59734 
1.61823 
1.63910 
1.66018 
1.68130 
1.70249 
1.72369 
1.74566 
1.76236 
1.78501 


1.80913 
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Table 37* 


POTENTIAL ENERGY FOR THE 02 Ber STATE 


v V(cm™!) 
0 353.45 
1 1041.45 
2 1706.55 
3 2347.95 
2965.55 
3557.35 
4118.45 
4651.65 
5153.05 
5618.55 
6047.45 
6435.65 
6781.05 
7081.75 
7336.75 
7545.15 
7714.11 
7847.02 
7949.15 
8024.79 


8077.44 


*From TURNGPT 


Fin 
1.52754 
1, 46883 
1.45786 
1.43628 
1.41890 
1.40495 
1.39306 
1.37998 
1.36959 
1.36125 
1.35429 
1.34833 
1.34334 
1.33902 
1.33507 
1.33088 
1.32861 
1.32627 
1.32496 
1.32740 


1.33034 


r 
max 


1.68252 
1.75912 
1.81448 
1. 86469 
1.91241 
1.96052 
2.00873 
2.05627 
2.01829 
2.16443 
2.22607 
2.29561 
2.37479 
2.46677 
2.57613 
2. 70806 
2. 86245 
3.04965 
3.27710 
3.56375 


3.94952 
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MOLECULAR DATA NEEDED FOR TURNGPT (02 x*r7) 


x 


0 


i 


Table 38 


Gv(em=!) 


789.81 
2352.20 
3891.68 
3408.25 
6901.91 
8372.64 
9820.43 

11245.26 
12647,12 
14026 .00 
15381. 88 
16714. 79 
18024.69 
19311.59 
20575.48 
21816. 33 
23034.10 
24228.74 
25400.22 
26548.54 
27673.62 
28775.15 
29852 .56 


30905 .29 
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Bv(cm~‘) 


1.4375 
1.4215 
1.406 

1.3909 
1.3755 
1.3%01 
1.343 

1.3296 
1.3145 
1.298 

1.2833 
1.2678 
1.2516 
1.236 

1.2204 
1.2046 
1.1870 
1.1711 
1.1550 
1.1390 
1.1226 
1.640 

1.0904 


1.0744 
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Table 39 


MOLECULAR DATA NEEDED FOR TURNGPT (0, Br, ) 


v Gv(cm™!) Bv(cm™!) 
0 353.45 0.813 
1 1041.45 0.798 
2 1706.55 0.785 
3 2347.95 0.770 
4 2965.55 0.754 
5 3557.35 0.735 
6 4118.45 0.719 
7 4652.65 0.703 
8 5153.05 0.680 
9 5618.55 0.655 
10 6047.45 0.627 
. 11 6435.65 0.596 
12 6781.05 0.562 
13 7081.75 0.525 
14 7336.75 0.485 
15 7545.15 0.443 
16 7714.11 0.398 
17 7847.02 0.354 
18 7949.15 0.306 
19 8024.79 0.258 
20. 8077.44 0.208 
21 8109.38 0.159 
: 22 8126.28 0.109 
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Table 39 (cont'd) 
CONTINUUM VALUES 
9150.0 1.31 
9560.0 1.305 
9970.0 1.300 
10380 .0 1.295 
40780.0 1.290 
11200.0 1.285 


11640.0 1.28 


Table 40 


POTENTIAL ENERGY INPUT DATA USED FOR CALCULATING THE FRANCK-CONDON 
FACTORS FOR O27 SCHUMANN-RUNGE BAND SYSTEMS 


V(cm7!) Thin rmax 
353.45 1.52754 1.68252 


Ie 


oO 


| 1041.45 1. 46883 1.75912 
2 1706.55 1.4578 1.81448 
3 2347.95 1.43628 1, 86:69 
4 2965.55 1.41890 1.91241 
5 3557.35 1.40495 1.96052 


6 4118.45 1.39306 2.00873 


8 5153.05 1. 36959 2.10829 


9 5618.55 1, 36125 2.16443 
10 6047.45 1.35429 2.22607 
1l 6435.65 1, 34833 2.29561 
12 6781.05 1.34334 2.37479 
13 7081.75 1. 33902 2.46677 
14 7336.75 1.33507 2.57613 
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7 4652.65 1.37998 2.05627 
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Table 40 (cont'd) 


v V(em™!) ‘nin "max 
15 7545.15 --- 2.70806 
16 7714.11 --- 2.86245 
17 7847.02 wee 3.04965 
18 7949.15 --- 3.27710 
19 8024.79 --- 3.56375 
20 8077.44 --- 3.94952 


7. PHOTOIONIZATION FRANCK-CONDON FACTORS 
a. No 
The four photoionization transitions considered were 


(1) No x1zy +Nt xr : 
(2) No xine +N an, 


ly? + .o.+ 
(3) My Se ee » 


let + .2.+ 


All five potentials used were generated using program TURNGPT and with the 
exception of the B state were generated from spectroscopic constants given by 
Gilmore (Ref. 21). The B state was generated by use of constants from Herzberg 
(Ref. 29). 

Table 46 contains the spectroscopic data used. Tables 41 through 45 


contain the generated potentials. 
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Table 41* 


POTENTIAL ENERGY FOR THE N> xr STATE 


vy V¢em™1) Fin *max 

0 1175.49 1.05485 1.14544 
1 3505.14 1.02644 1.18417 
2 5806.91 1.06820 1.21288 
3 8078.91 0.99410 1.23755 
4 10322. 88 0.98242 1.25993 


5 12538.14 0.97236 1.28081 
6 1.4724.61 0.96349 1.30064 
7 16882 .22 0.95552 1.31968 
8 19010.90 0.94829 1.33813 
9 21110.57 0.94164 1.35612 
10 23181.15 0.93550 1.37375 
11 25222.58 0.92978 1.39111 
12 27234.78 0.92443 1.40825 
13 39317.67 0.91940 1.42522 
14 31171.18 0.91464 1.44208 
15 33095.24 0.91014 1.45884 
16 34989.77 0.905 86 1.47556 
17 36854 .69 0.90178 1.49224 
18 386 89 .94 0.89788 1.50893 
19 40495.43 0.89415 1.152564 


20 52271.10 0.89056 1.54240 


*From TURNGPT using Gilmore's constants. 
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Table 42* 


POTENTIAL ENERGY FOR THE Nt xr STATE 


I< 


20 


V(cem7!) 


1099.54 
3274.34 
5416.60 
7526.22 
9603.09 
11647.12 
13658.25 
156 36.45 
17581. 70 
19494.00 
21373. 38 
23219. 89 
25033.61 
26814.62 
28563.05 
30279 .03 
31962. 72 
33614.31 
35233.99 
36 822.00 
38378.58 


Tain 


1.07397 
1.04461 
1.02575 
1.01120 
0.99917 
0.98884 
0.97978 
0.97170 
0.96440 
0.95777 
0.95169 
0.94611 
0.94096 
0.93619 
0.93178 
0.92768 
0.92389 
0.92037 
0.91712 
0.91412 


0.91136 
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r 
max 


1.16768 
1.20792 
1.23791 
1.26385 
1.28752 
1. 309 76 
1.33100 
1.35156 
1.37161 
1.39130 
1.41075 
1.43002 
1.44919 
1.46831 
1.48743 
1.50659 
1.52582 
1.54515 
1.56462 
1.58424 


1.60404 


*From TURNGPT using Gilmore's constants. 
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Table 43* 
POTENIIAL ENERGY FOR THE No A’n | STATE 

v V(cm7!) “min Tax View) 
0 947.96 1.02656 1.22743 0.1175 
1 2821.38 1.09570 1.27158 0.3498 
2 4664.80 1.07615 1.30488 0.5783 
3 6478.22 1.06118 1.33368 0.8031 
4 8271.64 1.04887 1.38479 1.0242 
5 10015.06 1.03838 1.38479 | 1.2416 
6 11738.48 1.02917 1.40847 1.4553 
7 13431.90 1.02097 1.43137 1.6652 
8 15095. 32 1.01357 1.45371 1.8714 
9 16728. 74 1.00681 1.47563 2.0739 
10 18332 .16 1.00061 1.49726 2.2727 
11 19905.58 0.99487 1.51867 2.4678 . 
12 21449.00 0.98954 1.53995 2.6791 
13 22962.42 0.98455 1.56116 2.8468 
14 24445. 84 0.97987 1.58235 3.0307 
15 25897.26 0.97546 1.60356 3.2109 
16 27322 .68 0.97130 1.62484 3. 3880 
17 28716.10 0.967356 1.64622 3.5601 
18 30079.51 0.95361 1.66774 3.7291 
19 31412.94 0.96003 1.68943 3.8944 
20 32716. 36 0.95662 1.71132 4.0560 


*From TURNGPT using Gilmore's constants. 
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I< 


V(cem-!) 
1197.15 
3588.66 
5887.46 
8147, 86 

10344.26 

12467.06 

14508.06 

16469.16 

18291. 36 

20018. 26 

21614.96 

23094. 8 

24466 . 26 

25742.56 

26938. 86 

28065. 46 

29132.56 

30148.06 

31114.06 

320 36 .06 

32858.06 


Table 44 


POTENTIAL ENERGY FOR THE N> 


Tain 


1.03331 
1.00546 
0.98795 
0.97409 
0.96323 
0.95501 
0.94720 
0.94043 
0.93404 
0.92856 
0.92315 
0.91867 
0.91496 
0.91145 
0.90846 
0.90544 
0.90303 
0.90104 
0.89893 
0.89687 
0.89012 
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+ 
Br STATE 


r 
max 


1.12304 
1.16196 
1.19199 
1.21811 
1.24330 
1.26903 
1.27401 
1.32000 
1.34765 
1.37643 
1.40798 
1.44118 
1.117632 
1.51208 
1.54832 
1.58452 
1.62095 
1.65763 
1.69439 
1.73160 


1.78474 


"(ev) 

0.1484 
0.4424 
0.7299 
1.0101 


1.2824 


1.5456 


~ 1.7986 


2.0405 
2.2677 
2.4818 
2.6797 
2. 8632 
3.0332 
3.1914 
3.3397 
3. 794 
3.6117 
3.7376 
3.8574 
3.9717 


4.07% 
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Table 45* 
POTENTIAL ENERGY FOR THE Nt crt STATE 

y Voem@!) Fin Fmax (ev) 
0 1034.08 1.21465 1.31124 0.1282 
1 3085.06 1.18158 1.34967 0.3824 
2 5111.94 1.15962 1.37772 0.6337 
3 7113.28 1.14238 1.49182 0.8818 
4 9087.64 1.12798 1.42383 1.1266 
5 11033.58 1.11558 1.44462 . 1.3679 
6 12949.66 1.10471 1.46470 1.6054 
( 7 14832.44 1.09508 1.48437 1.891 


166 86.48 1.08646 1.50386 2.0687 
18504, 34 1.07874 1.52333 2.2941 
10 20286 .58 1.07182 1.54294 2.5150 


eee eee 


*From TURNGPT using Gilmore's constants. 


b. 0> 


The photoionization transitions considered were 


(1) 02 Xx3zF + oF 


(2) 


(3) 


(4) 
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The seven potentials used were generated by program TURNGPT. The 0, XE state 
was generated by using data from Wallace (Ref. 24). This is the same potential 
that was used as the lower state of the 0, S-R system. The six states of ot 
were all generated using the spectroscopic constants given by Gilmore (Ref. 21). 
These data are listed in table 53. The potentials are listed in tables 47 


through 52. 


Table 46 


SPECTROSCOPIC CONSTANTS FOR No PHOTOIONIZATION STATES 


We  Wexe «= eye = Weze =e ae ye re 
Ny X!x* state 
2358.07 14.188  -1.2472 0.0 1.9987 0.0178 0.0 1.0976 
Ne x2x* state 
2207.17 16.146 -2.85-2 9,2-" 1.9258 0.01743 -1.6%7" 1,118 
Ne Ar state 
1903.42 15.0 0.0 0.0 1.748 0.02 0.0 1.117 
Ny Br? state 
2418.7 22.53 ~6.7-1 4,072 2,083 0.0183 -1.6573 1.075 
Ny cert state 
2073.7 10.97 -2,47! 0,0 1.5124 0.0017 -7.5°" 1,262 
Table 47* 
POTENTIAL ENERGY FOR THE 05 X?n STATE 
y V¢em"!) Fin Fnax (ev) 
0 934.12 1.07839 1.17354 0.1158 
1 2777.56 1.04947 1.21552 0.3444 
2 2587.94 1.03116 1.247.4 06.5678 
3 6365.26 1.01714 1.26466 0.7891 
4 8109.52 1.00561 1.29990 1.0054 
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Table 47* (ccat'd) 


y Vien) Fmin Fux (ev) 
5 9820.72 0.99575 1.32370 1.2175 
6 11498, 86 0.98710 1.34652 1.4256 
7 13143.94 0.97939 1.36865 1.6295 
8 14755 .96 0.97241 1.39030 1.8294 
9 16334.92 0.96603 1.41161 2.0252 
10 17880. 82 0.96016 1.43268 + .2168 
ll 19393.66 0.95471 1.45362 2,4044 
12 20873.44 0.9962 1.47448 2.5878 
13 2320.16 0.94484 1.49534 2.7672 
14 23733.82 0.94034 1.51624 2.9424 


*From TURNGPT using Gilmore's constants. 


Table 48* 


POTENTIAL ENERGY FOR THE Of Ay | STATE 


v Viem™!)  Fnin Fnax (ev) 
0 446.08 1.34603 1.48386 0.0553 
1 1317.68 1, 30607 1.54794 0.1634 
2 2161.88 1.28114 1.59762 0.2680 
3 2978.68 1.26225 1.64192 0.3695 
4 3768.08 1.24682 1.68348 0.4672 
5 4530.08 1.23364 1.72353 0.5616 
6 5246.68 1.22209 ', 76277 0.6505 
7 5971.88 1.21173 s.80167 0.7404 
8 6651.88 1.2024” 1, 84057 0.8246 
9 7304.08 1.19358 1.87975 0.9055 
10 7929.08 1.08541 1.91945 0.9830 


*From TURNCPT using Gilmore's constants. 
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Table 49* 
POTENTIAL ENERGY FOR THE Of b*n STATE 
V(cm7!) Thin Tmax View) 
594.11 1.22498 1.34439 0.0736 
1756.70 1.19630 1.39967 0.2178 
2885.11 1.16869 1.44240 0.3577 
3979.34 1.15234 1.48038 0.4933 
509.39 1.13900 1.51592 0.6248 
065.26 1.12765 1.55006 . 0.7520 
7056.95 1.11772 1.58341 0.8749 
8014.46 1.10886 1.61637 0.99% 
8937.79 1.10082 1.64922 1.1081 
9€26.94 1.09343 1.68218 1.2183 


10681.91 1.08656 1.71545 1.3243 


*From TURNGPT using Gilnore's constants. 


Table 50* 


POTENTIAL ENERGY FOR THE 0+ an STATE 


Vicm=!) min T max " (ev) 


515.25 1.32207 1.45020 0.0639 
1530.05 1.28390 1.50774 0.1897 
2425.29 1.25993 1.55142 0.31300 
3497.64 1.24174 1.58969 0.4336 
4450.21 1.22690 1.62500 0.5517 
5382.00 1.21429 1.65846 0.6672 
6293.01 1.20332 1.69071 0.7802 


7183.24 1.19358 1.72214 0.8906 
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Table 50* (cont'd) 


v Vom!) ‘min *max 

8 8052.60 1.18484 1.75304 
9 8901.36 1.17690 1.78360 
10 9729.25 1.16963 1.81397 


*From TURNGPT using Gilmore's constants. 


Table 51* 


POTENTIAL ENERGY FOR THE 0+ c'L” STATE 


y vier!) Fin *nax 
0 1195.50 1.14677 1.23086 
1 3559.50 1.12238 1.26895 
2 5887.50 1.10740 1.29783 
3 8079,50 1.09628 1.32307 
4 10435.50 1.08740 1.34627 
5 12655.50 1.08006 1. 368916 
6 14839.50 1.07383 1.28917 
7 16987.50 1.06847 1.40954 
8 19099.50 1.06382 1.42945 
9 21175.50 1.05974 1.44904 
10 23215.50 1.05916 1.46839 


*From TURNGPT using Gilmore's constants. 
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Yew) 
0.1482 
0.4413 
0.7299 
1.0141 
1.2938 
1.8390 
1.8398 
2.1060 
2.3679 
2.6253 


2.8782 
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Table 52* 
POTENTIAL ENERGY FOR THE 03 "ro STATE 
V(cm !) 
360.50 1.28272 1.43501 0.0447 


Thin Tmax . (ev) 


1066.50 1.23627 1.50491 0.1322 
1752.50 1.20657 1.55759 0.2173 
2418.50 1.18361 1.60411 0.2998 
3064.50 1.16440 1.64740 0.3799 
3690.50 1.14791 1.68880 0.4575 
4296.50 1.13310 1.72910 «0.5327 
4882.50 1.11970 1.76878 0.6053 
5448.50 1.10730 1.80821 0.6755 
5994.50 1.09572 1.84768 0.7432 
6520.50 1.08478 1.88740 0.8084 


*From TURNGPT using Gilmore's constants. 


Table 53 
SPECTROSCOPIC CONSTANTS FOR 02 PHOTOIONIZATION STATES 


Wexe Weye Weze Be ae re 


03 X2n_ state 
as: ae 


1.67220 0.01984 1.12270 
os a‘n state 
eS Shia 
1035.69 1.10461 0.01575 1.38126 
Os A2n state 
a u 
1.06170 0.01906 1.40890 
0; bYs7 state 
pera SaaS 


1196.77 1.28729 0.02206 1.27953 
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Table 53 (cont'd) 
_ We Wexe Weye Weze Be ae 


ee < 
0, 7Z" state 
ei 


726.00 10.00 0 0 1.15000 0.01700 


Of c#r7 state 
Se! Sore 


2400.00 18,00 0 0 1.50000 0.02000 
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ls 


0.0 


0.0 


[3 


estimated 


estimated 
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SECTION III 


TABLES OF FRANCK-CONDON FACTORS 


This section contains tables of the Franck-Condon factors calculated from 
the data given in Section II. In Appendix I is a listing of program TURNGPT 
which was used to generate the input data for program FRANKON. The Franck- 
Condon factors have been calculated by program FRANKON which is listed in 
Appendix II. 

Table 94 contains a list of R-K-R Franck-Condon factor tabulations known to 
the authors at the time this work was started. This list includes only transi- 
tions of interest for atmospheric diatomic molecules. 
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Table 94 


TABLE OF FRANCK-CONDON FACTORS BY VARIOUS INVESTIGATORS 


Molecule System Re ference i 
CN Red (A-X) 33 
CN Violet (B-X) 33 
co 4(+) (A-X) 33 
co* Comet tail (A-X) 34 
co* (1- ) (B-X) 34 
No (1+ ) (B-A) 33, 35, 

; 36, 37, 

29 

No (2+ ) (CB) 33, 35, 
36, 37 

No Vegard Kaylan (A-X) ee 37, 
No y bands (B'-B) 3% 
No y bands (a-a‘) % 
No y bands (w-a) 36 
No LBH (a-x) 37 
No Tanaka (c-x) 37 
No B-H (b', b-x) 33 
NS Meinel (A-X) 33 
Ne (1- ) (B-X) 33 
ND (2- ) (CX) 33 
NO 8 (A-X) 33, 38 
NO y (B-X) 33, 38 
No § (C-X) 33 
NO € (D-X) 33 
NO B* (B'-X) 33 
NO y' (E-X) 33 
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Molecule 


NO 


NO 


NO 


NO 


Table 94 (cont'd) 
Syster 
Lagerquist-Miescher (G-X) 
Heath (C-A) 

Feast 1 (D-A) 
Feast 2 (E-A) 
Feast 2 (E-C) 
Feast-Heath (E-D) 


Schumann-Runge (B-X) 


Reference 


33 
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APPENDIX I 


COMPUTER FROGRAM TURNGPT 


A listing of the computer program used to generate the RK potentials at 
AFWL is included in this appendix. The input data needed is described and a 
sample data deck is also included. 


1. INPUT DATA 


The reduced mass of the molecule or the mass at each atom is needed. These 
can be based on either C)2 or 0)¢. A good source for reduced masses based on 
0;, is Herzberg (Ref. 31). 


The energy levels of the molecule can be ccastructed intemally in the 
program if the spectroscopic constants we, wexe, weye, weze, wete are read in. 
These are used in the equation. 


E = we(vts) - wexe(vtls)? + weye(vt}s) ° 
+ weze(vt;)" + wete(vts)° +... 


Herzberg (Ref. 31) has a large tabulation of these. Alternatively, one can use 
the G, obtained from spectroscopic data. Wallace (Refs. 24, 40) has an exten- 
sive tabulation of these. 


One also needs the relational constants, that is, either Bee ae» Soe Ye? 


€,» or the By for the various vibrational levels. Again Herzberg (Ref. 31) is 


a good source for Be» O13 Ye and Wallace (Refs. 24, 40) is a good source for 


e 
the BY. The final bit of data about the molecule that is needed is the equilib- 


rium internuclear separation. This can be found in Herzberg (Ref. 31). 
2. DATA DECK 


The data deck set-up procedure needs very little explanation since Zare 


has made extensive use of conment cards throughout his program. 
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Card Number 


1 


= 3 


6a 


6b 


ITEST 
11 


FMT 
8A10 


IIMS, ZMAS1, ZMAS2 


IF ,2E10.0 


IQHK, N 
214, 4x, 6A10 


WE, WEXE, WEYE, 
WEZE, WETE, 5E10.0 


G(I) 
Format in columns 
13-72 of Card 5 


IOPEV 
11 


IBHK, N 
214, 4x 


118 


1 if problem follows 


Title of the problem with carriage 
control in column 1] 


IIMS = 1 mass based on C = 12 
IIMS = 2 mass based on O = 16 


ZMAS1 mass of first atom (or 
reduced mass of molecule) 


ZMAS2 mass of second atom (or 
blank when using reduced 
mass) 


IQHK = 0 using a tabulated Gv 
curve. Format for Gv in 
columns 13-72. 


IQHK # 0 using constants to 
generate Gv curve. 
Hollerith text in columns 
13-72. 


N = number of vibrational levels 
to be used 


Use when IQHK #0. These are the 
spectroscopic constants that 
generate the energies of the 
vibrational levels. 


Use when IQHK = 0. These are the 
Gv of spectroscopy. 


IOPEV = 0 if zero point energy to 
be found by extrapolation 


IOPEV = 1 tabulated Gv curve to 
be used unchanged 


IOPEV = 2 if Gv curve is to be 
constructed from constants 


IBHK = 0 using tabulated Bv curve. 
Format for Bv in colums 
13-72. 


et ET 
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Card Number 


Card 8 (cont'd) 


9a BE, ALPHAE, 
GAMMAE, DELTAE, 
5E10.0 

9b B(v) 


Format in columns 
13-72 of Card 8 


10 TOPA, BEQUIL 
14, 3F10.0 


ll RE 
F10.0 


12 VSTART, VFIN, 
HDED 
3E10.0 


13 IOPFG 
Il 


IBHK ¢# 0 constants are used to 
generate Bv curve. Holler- 
ith text in columns 13-72. 


N - number at vibrational levels 
(same as above) 


Use when IBHK ¢ 0; these were 
used to generate the Bv 


Use when IBHK = 0 


IOPA = 0 value of BE to be found 
by extrapolation from bv 
data read in 


IOPA = 1 value of BE given by 
BEQUIL 


IOPA # 2 rotational constants are 
used (leave BEQUIL blank) 


BEQUIL - spectroscopic constant bk: 
Intemuclear equilibrium distance 
Turning points are calculated fro 
vibrational level VSTART to VFIN j 


to steps of HDED 


IOPFG = O intermediate printout 
skipped 


IOPFG = 1 intermediate steps in 
an evaluation of the Kle! 


integrals are to be print 
out 


SROGPAT THEMCOT (INST pQUTPITT gs TAPES ZT MPLIT STADE KEOUTPUT oF ILMPL eFUNCH A 1 
Le TADE Ze incHy A 2 
3 SET ONO FAR CLOTTING OOTEMTT ALS coe A 3 
CAR CEE Ra er A 4 
C SATN FIMES TUSRING OOTNTS FOP A KALECULE BY PYQRFRG=KLE INSOer ES wre A 4 
AEC CONF SS OR SVE PP IMBERS AST ANGE TOIOQVS, A 7 
S mNeeD SOFCTELES CHOICE MPETWEER MASS UNITS RASED ON C12 = 126 OR A a] 
3 Ate = ihe ee 
c THE FOLLOM TNC COMMENT CARDS PESCRISCE THE DREDARATIA" OF NATA CARDS A 14 
r A 11 
AIYVENETON VSTPOT (TEs VEINIT Io HOFUCT) ; A 5 
PIMPS TOW NEECPRE Oe VECIAF Ls WAIIOPA) ©» VERIAN © SCFM De MVIZED)»s F&F A 17 
LIRTAL EATS De CLPAC2]) «© CTENC3) © QIMSE 2) © ZIMS( ADs PISAIA IL» ZFIENI ADS A 13 
PECEL TC 27% De KOON(K Es VEPN(S)e OTERYT (B00 GI2NDDe PVE2ZNODs FMT 1A)» A 164 
TUIEP MY e OMIM E POUL © PYAKEP0U) © DUFM ETB)» OVEL JON s GEEPOM)s RYTC2R A 16 
S710 S027 Ve WHATS ARVELIANDS FVELLIAINIs PVED(ANDs ALINIs VEIADs X¥KIH A IA 
i Bye YYCADe FHEC ADS MANADO XFL CL) 0 XFECTI® FLO 1)» C302)0 GBI3I0)s YY A 17 
AVv(8)e SICHOL) AP 
y COMMON WL eV TD oXm eo Vor ei TAN a TA ote KMIN GKVAY GINEC oMREG NISSEN oSotlfe A 19 
UES pMAXTT OF SCM eg 2M aE OWE ONE SFY E hE ZF gp WE TE SPE AL OHAF gCANMAFE QS HEL TA A 2a 
ZEA er aM eS THE ett aP ol o1CK otte 1 TZ HIMES oF OSL NE oKMVeAFQUTL e7OTEN oMLV oL A 21 ° 
MV ePFeVSTART VEIN eHNED A 22 
POSE SLeSP eS te Jha JHAE? A 24 
(a A 24 
.S FIRST CAPD Th OATS Ae A ONF IN COLIN 1 TF A POORLE™ FALL W'S. A 26 
€ VERY LAST CBRPR IN SATA MIIET OF A ALANK CARN eecee A 26 ty 
C a 29 
i COM TITUS A 28 
QFAN (8,8) PRCT A 29 
TF CLTETT) Fesre? A 34 
¢ A 4} 
£ NEXT T''O CAPHE IN DATA SAVE MATTE AF PROBLEY IN COLUMNS LekNs A 632 
i< WHEQE CAPRIACE CONTOCL TS IN COLUM Te NEXT CARD HAS TEMES AND A 28 
G MASSES CF THE THOM ATOMS. OF TIVES AND SENUCFA MASS ITN FIRST MASE A 4 
kK FIELG WITH SECOND FLANK. A 48 
¢ TIME w de VISS UNITS SOF RACED ON C12 = 176 A 34 
€ TIMMS 8 Pe VESS IMIITS ASF RESFO OF ATE = The A 37 
c A aq 
é DEAN (HOEL)D CFTC Tel atelA) A 39 
OTE~ (ROEDY TTC THAT IID m4 
A 4) 
é MPYT CME Tt METS FOE POQuitg ANI SHIMRER QE LEVELE, FAT A 6? 
C I! Te ECR ET. AS WELL AS CHE FORMAT STATIVENT vonleH CANTROLS THE A 4&7 
¢ RESINS SE THR LEVELS (3% COLUMN PAR TP me AD EVAVALEe (4F 1&6 eM de A 44 
iG [Wee 2 S TF THE FOLLOWING CARDS CONTATAN THE TaAntiLatre) % CHovyEr, A 45 
G TF A PAR TICULA® VALUE TS TO BF INTEFPPOLATEXs SET GIT) = -1Me A 464 
C ff 1° 1& SFSTRED TN USE CONSTANTS TO GENERATE THE FNTIPE A 47 
C G CURVE Tate 2ST HOT FF £QUAL TO ZERO™ A 64a 
S * GALLI TITH TRYT SHERT RE DIEM CHED TN COOLING LAMM 4 49 
AT FT WELL PE PRIPT EM» Fetes Che: oe 
rc THe Mervyy Cor fe Geet £ FACE CONTA ge 5 5} 
HE gr Myre’ FYE SKE PE CAMs WETE, A e&9 
- s 89 
DEAT CR QRT) PQ Me (NTEOVT( Ty elelot) sos 
| rroNe? 4 88 


120 


an DNNADAANIFNNNADANNNANANNS w 


NANANANNNANDANNNNDANNAN:S 


NAY 


AN 


IF (ICHK) 2049? 

READ (5054) WE oWEXF gvEVE gE 7F gif TE 
GO TO § 

READ (SeOTFRMT) (GICT) oTzleN) 


THE ZERO POINT FNEFGY “AY PE FOUND BY FYTRAPOLATION, OF IT “AY FF 
FIXED BY THE SER FROM HIS TAPULATED GIV) CURVE. 

NEKT CAPD CONTAINS INPEVe 

SET IOPEY = Se IF THF 7FRO POINT FMERGY IS TO PE FOUND Ay 
EXTRAPOLATION» AND THE GV) CURVE ITS TA FE BAISEN 02 LOWERED 
ACCIONINGL Ye 

THIS TS PARTICULARLY USEFUL IF DELTA GV) DATA TS TO FE USED. 

SET IOPFV = le IF THE TARULATFD GtV) CURVE TS TO RF USED UNCHANGED 
SET IOPEV = 25 IF THE G(V) CURVE IS TO BE CONSTRUCTFD FROM CONSTAS 


QFAN (5e89) TOPEV 


THE NEXT SFT OF CARDS CONTAINS THE PV CURVE WITH THE SAME 
RESTRICTIONS AS FOR G CURVE ASOVF. NUMBFR OF VALUES HERE MUST AE 
FQUAL TO THF “SUMRER Aanyver 


RFAD (59572) TRUK Ne (MUEMT(T) e721 96) 

TIEARe2 

TF CTBHK) Be 70h 

READ (5554) PF eALPHEE gGAMMAF DEL TAF oF OSLNE 
50 TO R 

OFAN (5eMUFYT) CRVT CT el=leN) 


THE ROTATICNAL CCNSTANT PF MAY FE FOUND BY FXTRAPOLATIONs OR IT 
vhY AF -TYFS OV THE HSER, 

NEXT CARD CONTAINS TOPA APD PEQUILe 

SET ICPA = Cy IF THE VALUF OF PE IS TO BE FOUND PY EXTRAPOLATION 


FROM THE RV DNATA READ INeo LEAVE RE QUIL BLANKe 

SET IOPA = le IF THE VALUF OF PE IS TO RE GIVEN BY PEQUILe 

SET TOPA & 2e TF RCIAITONAL CONSHANES ARE SED, LEAVF PFQUIL 
AL ANK » 


QFAM (5082) TODA SPEQUIL 


EAD IN MIGHF2 TERNS If EFEYEANSION MVeHVeCOHIOLING CONSTANT Ae LANDL, 
GAMMA OFA MULLIKFN FORMULA NUMPER CF J LEVFLE TO BF CALTILATED Noe 
AND IND (21 2 72) TO INDICATE FXPANSTON FOR’ 1 FORD SINGLET SIGMA 

9 FOR NOUPLFT EL 172 3 FOR NONRLET OF 272 

INN25 FL NOVELFT SIGMA 

IND F2 MOVARLET ST GPA 


REAN (55S) HVeLAMeNJeGAMe IAD 
PFAN (8 e546) (PVIT)eTeleN) 
READ §Be (ACT Peeled oN) 

MEAN PATATIOVAL LEVELS J 

PFAM (8&7) CUOT efeleNS) 


MEXT DATA CAFM CONTAINS REe THE FMTESNUICLEAD SSUTLISRTUM DISTANCT, 


12] 


~~ D> ow: 


i a ae 
st > 


aT. 


nn 


NIAINNAN 


NAANAANNN 


NNNAN 


MEAR (B,Fe) OF 


WEXT DSTA CAPD CONTAING VSTAPTOVEINSAND HNER WHERE 
TURNING POINTS APE CALCULATED FROM VIBRATIONAL LEVE'. VSTART TO 
VEIN IN STEPS AF RDFA Fete FROM O08 TO 2008 IM STEPS OF 1005 


BEAN (G,87) YETART SVE INGHNEN 


WENT CAS" CANTO INS yaAcer, 

SET TOPE = Ye TF INTFRIVENTATFE STEPS IN THE EVALIIATICN OF THE 
KLEIN ACTEM I" TEGPALS ARF TO °F PRINTEDs 

GET TOPOFG 2 A, TF INTERMEOTATFE PRIMNT<QUT 16 TO AF SKITRPPFN, 
NORMALLY THF '!ISFR SHALILA SET I99FG = De : 


MEAN (8 —h%) TNPEG 

PFAD £0, ©? 

RFAD Sle (270Tof2107) 

THIS TERY INATES COMMENT CARDS CN THE PREFOARATION OF DATA» 

TF ND FURTHER PRORLF'S FOLLOWs 2PEMEYBES TO ADD A RLANK CARD TO 
THE NATA NECK, 


ATRA(ZY2EHANCST, 
QTENC2»HAHTIC 

AIVS 01) e6HlC1 2212 
AIMS (02) s6HO1As 14 


POINT MEANTYNC 


WRITE (HobTV COUTCPPeleteté) 
WRITE (HeS1) CEMT( LT) of eletA) 


PRINT THE MACEES ANN THETO HINITS, 


TF (7MASPDY 2% o1%eN 

ARETE (KoA 7) DEMSELT SD e2MAL] o2AASP 
GOTO. VI 

WRITE (Ao&2) TIMECTINGS) oe2VAS) 


SINT GOV PATLe CR VIORATICMAL CONSTANTS 


TE (LQeMCy TPelrel? 

WRITE (Ke TFONT) 

WRITE (he64) WE etEXF oWEVE gWwE7E WETE 
G0 TA 34 

WRITE (HoeAEd STENCT TEN) 

WRITE (6 056) O71 OT oe Te) 0°!) 

WINTTE (A eKT) 


ORPINT WW TATA, OP RNTATICMAL C7 NS TANTS 


TE CITY Peete gtE 
WRITE [Re OEY) 
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A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
& 
A 
A 
A 
A 
4 
A 
A 
A 
A 
A 
A 
4 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 


18 


19 
29 


21 


NAAN 


NAIANN 


WRITE (6068) RE eALPHAFE pGAMMAE SDFL TAF SERSLANE 
G60 TO 17 

WRIT’ (6969) NTEN( TER) 

WRITE (607%) CPVICT) stellen) 

WRITE (K967) 


PRINT HEARING 
WRITE (A081) CEYTODpele1014) 


TF (7°AS2) 1P 91519 

PMipe 7PPAS] 

60 TO 2¢ 
PMUEZMASTEPYAE IS (QMASTEZ2NAT2) 
ZIMS CT) £1 6% 

7IMS (02) 2—e999K 794 
ZTRACL I E10% 

7TRAC2)21 eBROTES 
2MUSZMUFZIMS CT INS) 

STENT LP E7MUBAALIAAOE2 
ZIFNOC2VE7UH1eABLOS26F M7 
PTEN( 3) e721 FIVTITAFP 
FACMEE0 e201 709/2MU 
X1(1)2028 

NO 21 Teen 

tleI-1 

VICTPSYTCTL 412% 
CONTINUF 

7EFOF0,% 


CALCULATE THE Av TO BE USED 1% LATER COMPUTATION CF PCTENTTALS 


J1206 

AO 4@ LJueleNJ 

JG2I(L J) 

LAM2=LAMFL AY’ 

TNMh#& FOR TRIPLET Sparta STATE 

BO TO C722 oP ae Pha AMe AWA gA*MeMM—e IMO YQy YON 


SINGLFT SIGMA "LFECTRANIC STATES 


J? 2 jIG*(J44)) é 

AO 2% Tele 

OV ET PEQVECT a2 eMrVIP IF IISUVE IDO SD 

PVE SAVIIT IPE S2—OV CLIP IDB ID SHV ID HIDE ID 
FIL VEGTIL per ver 

COUT ZGtT) 

CONTINUE 

GO TO 4] 


QOWRLET PT 1/77 FLFCTRONIC STATS - HILL “89 vAN VLECK= 


FI 

SHAPED = (JbagSyarr, 
Jte2# he} 
JleJGe (ad t8 OT 9417465 


>Pebpre > 


» 


re 


PP PP PP > PD fe 


rPPrPrPrrpPPrePrre rr rYrrpr> ppp >» 


p> th > 


165 
164 
167 
140 
149 
U4 
Le 
172 
17:9 
174 
175 
174 
177 
178 
179 
180 
181 
1A? 
1843 
184 
1A5 
186 
1A7 
yra 
129 
19Nn 
19} 
192 
}Q4 
194 
1986 
196 
197 
198 
199 
IAN 
ant 
209 
IAA 
ANN 
SAG 
onk 
ART 
2ne 
pao 
A1N 
11 
21? 
?}2 
P14 
C8 
214 
237 
21a 
219 


AQ PR Tale 
MARVICTISA TEAC TIFYISAT/OENEOV(T) 

AVITVERE (L(G et JHAFIOVIS (VI me) MLAMI) BERRA] / JI Het NE ILE H TQ = QSH 
sy 

AVE PREC MAF Del ANISH (Lg IHAF OY #1 Late) HLAN21 O98 aD UH HL +0586 
1EV8( Jo=05) 

AUT EGTELyenve 

rE(TpaGd) 

ecaurres tre 

80 7") 43 


NMOVRLET PL 37? ELECTRONIC STATE —- HILL ANA VAN VLECK oe 


2 

JHAE Da ( JL et peer, 

Jie JGelelL yee sineyy#eS 

J322* 144} 

AA OF Ts) eX 

MAOVITCTISALEACTISYI SALT /°ER 2OhVv(T) 

PVITIEPH LL e+ (Ge HIHAFIOVIF(V Labo) FLAM2) FH meh mle / JI (GENE ILEEIC +A5H 
1GA™) 

CVE eB (JHA MeL EOS RH ( gS IHAF DOV IF (YY mbe) FL AM2) OHS HN H( Ose) HEM, 
LeeS AGA ( JOO] ,%) 

GOL =ZGIeL peeve 

cle ee | 

CONTINUE 


GE TAN o4 1 
TRIPLET SIGMA FLECTFONIC STATE 


BalS COUPLING PAMAYFTED LAMA ASA 16 FMPTY FAD THIS CALCHILATION 
SI PRIAGPOA* eK (PUVEICAL ) 

TFAzy -9- FY] gfF M22? wee FP off PAI we FI 
pyrex 

POA TA (7% eA eI7IA~SS [FO 

| 

JAsI* Jae} 

'PeJant ner) 

Of we Le le® 

MPEEY I(T ISAT SAC TPENENV(T) 

JAG2( 27# J4eTyp Hn 

TODS ISTH ILM KP POAT ODBAT EP 
TEMSHMAY] O% O0TSR 

Tome ons Ss 

PV OT) = 22 RON IPS] SII (PHP RDHEILRERSTEQSGL") 
PVE SAY JIAA BIDE IPG ILMP MALH-TSNSC ANH ( S44) ) 
GOUT AGL CTyoeve 

G8(ty=eac)) 

COMT [haar 

oy IS 44 

cr? 

J?2JL¥(J447) 

aN 427 ey ot 

MED P(PPTA LSA PESOS WTy 

AV (Ty ePa te ABD 


A 
A 
A 
A 
A 
A 
A 
A 
4 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
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A 
A 
A 
A 
A 
A 
4 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
4 


AVYF=EPO# IJIN IPH ID 

ACT GIL peeve 

GPCi)=2G01) 

CONTINUE 

en TT. 41 

a] 

J22J4(544)) 

JIS7* J+] 

AO 26 TeleN 
MPSAVICT IEP Tes ( Tyre rviy) 

JGR2(2# JU—]) RO 
TSODBIG2H# SGM GAL IAL H9#A] SP 
TSMNEAMAYIT(% ee TS9) 

TSQ2( TSN) HAE 

AV IT) 2RH2 #8 I72 41/598 (eR OPH IGABASTSOMGAM) 
AVE BR #I2—J4R-AL4+TSO-GAV# JU =N8 2% 7 
GOD eG Cl +eve 

6B(1)26(1) 

CONTINUE 

GO TO 4} 


PrPp> PPpPr p> 


> 


NOURLFT SIGMA FLECTRONIC STATE 


AGA" FO F&F THIS CALC ILATION 

J US ACTUALLY K SINCE “VE HAVE CASE © COUPLING 
Fl 

JYs2t Jae} 

JP 254% (544+1) 

JleJ4 

NO 36 Tele 

PRaAVICTIS CACTI ENS°V(1) 

AVP EPH (1 ere PIN IP/ARKC/ (7g HIIHA)) 
PVE SR IP eDH IPH IPRKRCH III, 

G(DeGiel peeve 

GP(T)2GC1) 

CONTINUE 

40 TO 4) 

F? 

J722*IJ44) 

Je2Ja*(J&e)) 

J} 2Jh 

AX 120 T=) gf! 
MeQVTCTIICHACT EVE (T) 

AVIT PERE 1 ema HDF I2/P40/(74%44%I52)) 
AVY FSERH IPSC H ILI eR E IPH IP 
GUTeCTCTysenve 

G2 T)=sG0T) 

CONT DUE 

60 19 8} 


A 
A 
A 
A 
A 
A 
A 
a 
4 
A 
A 
f 
A 
A 
ra 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
a 


> 


?P PPP Pe > > P te 


te 


NOUALET PL FLECTISAMIC STATF = ALMY SNR HODGFALL = 


AH VOPT If ThMVFRTFN FPe(aseleK% 4) /7 AND © la(4)-seee]/* 
Kk 1S RFAD IN MiIcTt COMVERT TO J IND27 FOR F2 IN 29 FRO 7) 
JEPRIGRAM EK CPIIVSICAL) JSTCOROGRA* PS IIPHYST CAL) 


18 


rpprpPp > 


ao 


> 


IVNANNEF NG 
- 


ame) 


42 


ON NEN 


JT et(Libeleal pet TNnhe1)#.8 

J2e Ji #145} 

SO 4D TeleN 

ABRAVICT) 

BeNvV(1) 

CaAl(t) 

FaLAvVNAsévsry 

reC/a 

Fat ye hy say 

F25/8 

A] aha Pane JI an 

QARAL PHS wae 

DP aLOF R(GP IP HLH DEE SESE ( DH IIR BS LH ( 2H SIMA K OF RFE 
AA wl 25 *F ROME 41 ol S29 78) OAD) ERE 
DATTIAL DERIVATIVE CF ALPHA = SP 
RPweS/PAH(RI4( I2—e TES H(GHF SFU He ( OH J2—) #7) ) 
TeIND FOR FL gPesINh FO FD 

OVI LT pth lola] peeINSer apn 

FE ae mM eM) ENO ( SPH IPH Qh HIPG1 IS 16) SPRAAS(@} DP SAITAND 

COLI SGT(T oer 

PACT) ently 

CONTINUE 

PRINT K AND J THE ROTATIONAL CUANTUM NUMBERSS FOP THIS PRORLEY 

ORINT Fle J(LJ) oJ1 


ORINT SV AY JST CALCULATEA 


HQTTS (Ke?) ATEMOT TSE) 
WRITE (Ho77) (OCT Vel zl oN) 
MPITE (Fe7K) QTFNOTIF®) 
WRITE (407k) (OVO To Tet ofl) 
PRINT Fle SILI) oJ 


EIN ZEPO PCTNT VAPTAPLES 


MINTS] 

BO 4? [eleAé 

yx (Tt) =7 

v¥ (TEP yy) 

vyvOLyss0T) 

CORT PAE 

eTGm(]payy (fh) 

VAQOER 21 

LU Seas 

CALL LECPRT (X¥eYYVONINT oP OPPER OPTS SIGH ANS) 
REGQUILEANSCTISANE (27) HOVE VOANS (39 4( 6275) 

CALL LEGPPT CYVoYVV OMINTONORNER oMOTS SICH OANS ) 
PPTEMZANS( LIANE (PVC SSIOENS (394 (075) 

AF Sag 1LOSARICAHS(REQUTLS7YAS) ) HME Sh 

CALL TXT COMTN gS MAY OUSTART AVE To HOFD ONFINGEV] sUeT OHM oT PHK SLOPE VIA 
LO4 RET 'TL ee IOPFS) 

ARAET KN ALNCcY oF isloe ain 4 Oy 

REQUIEES ANTEFIINAL MATAR CARNES AT ENN OF HATA DECK 
SISOCTATIA™ Emrocy - CAIRN 1 F°4? AD 

HISPHC ATTA PEasic T< -~ CANE ? Oia 93 


126 
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230 
331 
432 
332 


3346 


346 
336 
237 
439 
139 
449 
341 
162 
343 
344 
345 
346 
347 
34R 
440 
35” 
351 
352 
282 
356 
355 
456 
ART 
35" 
ako 
46n 
461 
162 
1643 
264 
qAS 
IAG 
247 
368 
24° 
370 
a7) 
3472 
ayva 
374 
375 
3476 
377 
a7_A 
379 
3ANn 
an} 
499 
any 
284 


43 


44 


45 


PRINT Sloe (FMT CT) eT=1016) 

PRINT P60 (NIT ol=el.7) 

APF aN7/2Cb66e%? 

PRINT 77e 7¥AS1T IZ NIE 

PRINT FAe JCLJ) oJ} 

PRINT 79 

PRINT 2&9 

7E29-% 

PRINT Rilo 7F e727 oPF o7~ oPFQUIL 
ASTAR=]H* 

PUNCH 5le (FMI CTP oT=1 08) 

PUNCH R22 

Lest 

NG 44 TeleNFIN 

G2eu(1) 

612G62/9966eC2 

63262/1997374039 

RMTeOMTe CT) 

RMABOENMAX( 1) 

QMTLERV1/6£&291 772 

RPMAL#RMA/ 2529177 

VLEV2YI(T) 

CHEVLFV-INT(VLEV) 

TE (CHeGT eee etNeCHelTee8S) CO TO 43 
PRINT R60 VLEVeGloRMI oe AeG2eASTAR 
GO TO 44 

L2t¢] 

tI=l-1 

ALeAvil) 

AlsA(L) 

DlePvtL) 

ALTEAVI(L) 

G2TeGI(L) 

FeGR(L) 

FIlsF/8946—6%? 

PRINT B59 VLEV eG oR OMA GP ePl OF] oF oASTAVeA2L o2TLeMloAloll 
PUNCH 836 [J] 9G? 

CONTINUE 

PUNCH 845 PF ZF 

00 45 Tel] oNFIN 

PUNCH Pe REIN(G TP eUCT) eRVAX (CT) oth CT) 
CONTINUE 

TF (LIeSTel) GO TO 4& 

OMe ep (NE TM pee? 

MM EN ZANAX1 (199 OM) 
XMAXERMAX(MFIND O16 

YYAKEN2E200OC, 

CALL MAPG EXMIN oe XMAXK oCeeYVMAK eo oleleselole) 
CALL ARSLINE (%eeMeelCevle) 

CALL ARSVECT (lesle) 

CALL ARSVECT (J e0%e} 

CALL ASVECT (%e0%e) 

CALL ARSPFAN (,24,%F) 

CALL CHAROPT (%9%0%0%o%} 

CALL SYMPOL C1IOHRADIUSCANGST2OMS E60) 
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46 


47 


4a 


CALL ARSPFAS (626 00%?) 

CALL SYMPOL (28HRKR POTENTIAL GEPERATEDS$.e) 
CALL ARSREAY ( 668 00%?) 

CALL SY™RIOL CFT) 

CALL PREOEAM (rhea) 

CALL CHAROPT (Co%o%el 0%) 

CALL SYMBOL CLOHFNERGV(ICMOeL Se) 
v01p29%1) 

VOPO 1D ERE 

VO2N0] PERE 

ND 4? [eleNFIS 

Kale} 

veiw pest) 45801) 
VELOT OK) sPYMINGG) 
VEZOCN OK) FAV AV OT) 

NPTS#MFIhe} 

CALL LINES (VEIL eVIT)ONPTS) 
CALL LINES (V02971) 0V61) oNPTS) 
NL aJILJ) 

CALL CHAROPT (%o%o% oD e%) 

CALL MIMRED (bg THAME 612) 
CONTINUE 

CALL FOALNE 

60 TA 1 

CALL EXIT 


FORMAT (11) 

FORMAT (PA1") 

FORMAT (1402?F17%6%) 

FORMAT (2 1604% 96A1%) 

FOR™WAT (SEl oe) 

FORMAT (EL %eGoTlolte lt e% ell) 

FORMAT (8F 1%) 

FORMAT (10173) 

FORMAT (F17—%) 

FORMAT (FIC —%) 

FORMAT (F997) 

FORMAT (2A1°) 

FORMAT (///6%4 THE “ASSES CF THE THO ATOMS» BASED CF Abeble ARE & 
Lp eH eH AND FIL eAS/SS/) 

FORMAT (/ 7/64 THE PERICES AES OF THE THN ATOMS, PASEN OM AGeSHe 
1 1S FIDeAs//) 

FORMAT (63He Tut Gly) Creve 1° CONSTRUCTED FROY THE FOLLOWING Iu 
TT DATA. // TH “Fos EU ge Ae BHe ENF FIC ee RHy WEVE FING 39h, WEZE 
2elVVereMhte WETE FE10—3) 

FOVU.AT (32H THF INPUT G VALUPSe FMFERGY IN AGelRAHe ARE GIVEN SFLOW 
1e//) 

FORTIAT CI¥Y TTF LL 07) 

FOT'AT (27 Ue GFYVALIIFS SF -10, ARF TO PFE FOUND AV INTEQPOLATI Ole) 

FOSPAT (64H THE Evtiyy CURVE TS CONSTRUCTER FROM THE FOLLOWISG INF 
LIT P27 Ne S/T FE = JOE YL eh ole ALPHAF SEL] eel AH, GANYAS 20 teh, 
21740 CELTAE =F ll ehel “He FPTLNE =F ll es) 

FQORAT (23m THe [ThouT AY VALIIESs FNFPCY 1% AGel MM» ARF AIVEM FEL 
Tte//) 

FCRMAT (IXIOFL1 66) 
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440 
44) 
442 
444 
444 
445 
446 
447 
44% 
449 
4590 
45) 
45? 
454 
454 
4586 
456 
457 
45R 
459 
464 
46) 
462 
463 
444 
4658 
464 
GAT 
464 
469 
479 
471 
472 
473 
474 
478 
ats 
477 
47R 
479 
aan 
uF] 
482 
4A4 
RG 
485 
GRG 
4°7 
4RA 
GAO 
4or 
4%} 
al] 
| ae ec 
AL 


71 FORUNAT 
72 FORMAT 
1GY IN 
73 FORMAT 
14 CORMAT 
75 FORMAT 
76 FORMAT 
77 FORMAT 
19.7e6H 
78 FORMAT 
79 FORMAT 
60 FORMAT 


(GHIKs o13e4H Se oF he!) 

(61H THE FOLLOWING NUMPERS AME THE “ONIFIE™ © VALUFSe FNFP 
eAGBel 7TH APF GIVER REL Owe //) 

(1X1°F1123) 

(46h THE FOLLOWING ARF THF YONTIFIENR By» FRFORY IN Abe//) 
C2¥1CF11e6) 

CLHO o6X e2OHDISSOCTATION PRODUCTS ARE o2A90) 

CLHO bX eLIHREMIICED MASS ee FS eb eAIHDISSEOCCIATION FNEQGY eek 
TAC ROF MG OI FY) 

CIRC oS XK ol THROTATIOMAL LEVEL ohXo2tKee Se Tv a MHUe oF Eel //) 

(2B eISHCALCLLATES MATA eh ToL IHENOTOIMENTAL MATAS/) 

(GH Ve eS eGHGVIFV) 0 AX eo THRIIN( ADS 2X opTHAYVAK (A) oo AV oRHGVIT/SEM) 


1 e2X eAHBVI1/CY) 0 AX oSHF (EV) AK THE (10M) GX 1H # oO op AHGV(1/CM) 9 BY AHA 


2vil/cn 
al FORMAT 
82 FORMAT 
83 FORMAT 
As FORMAT 
65 FORMAT 

IXeoF7el 
86 FORMAT 

END 


Ye eMHOVIL/C%) 0 3¥ aBHAVIT/CM) ep AXe HV) 
C2K oF Ve le2XoFVebKobK oF Teohe TX oF Veter XK oF Beh) 
(11H KVeFTRIAL ) 
(T1%F1%.3) 
(6F 1208) 
(LH FG ele 2X oF Deb ePX oF Ted o2X oF Tee PK oF Fete AK oF Reber oFbe 92 
4X Al ody oF Bel 0 3X FR eho rX oF Fede IK oF Pore IX0l2) 
CLH Fool e2X oF Meh 02x oF Tel a 2V oF Feuer K oF Vet 3PVeAl ) 
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495 
496 
497 
49a 
499 
BAN 
&n) 
8N2 
£93 
ENG 
BA5 
var 
507 
5A 
509 
810 
511 
512 
513 
814 
515 
516- 
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CUBROITISE VOR CRMTM oPMAM OVSTAOVEIN GHIE To MOEVL Ue IOHK eo IPHK eo IOPFV el 
PPA OUIL o{ TPF) 

PER FINS TUTS ING DCTNTS FIN COV) AND AV RATA. OR FROM CONSTANT 
AIMETSEOP SEM CANT Dw EV EIA), AVIDPAE DS? RICA TEYPIZAINADS 11099" 
Lye PTE 2OT Pe RAM E DA Ve VYO2* Je FLAY ® ZOC(9D)o MICED > T2050 THE 
2H 10 L540) 

SEMENTION VEIT DO BEIM De MOLZVC De CEVIAAMDS OVE29M Ds WELANde VOI? 
vdeo £GC27 pe FR MELD AGAUSIS ede KGAUS(S e2Od 0 VESTA LP ® VEIN( Te H 


2mED ide “2 ¢ De FLEE Le PNT CFDs Y¥CB De KFE LCL De SIGH] De® YY IADe | 
WHET De THAT De HTMAM DS THEM De IMEBA) 
AIVEVeTA PVP (2% be oVA-2 i Lalas | 


COMMON FeTE He FV er eT eMST ene” T eDeXMINeXMAM SINSCOVOPGNUSFI ONT NI 
LenS eV AMT TOPE e2¥ Ue DF oe HE OW EME oWE VE oWEZE oWE TE eo RE GAL OHAF pGAMYAF ODFELT 
2AE oVNIN 2G? Ko SPER ot il eV e7 of CK oHoK oHOF So FPSLNE o XX gp P®RRAND g7PTEN g VL 
BWolVere 

GFUNCFE OX eC OOO COOCME TE UX DOVE ZE POX DOME YE) OX DOWEXE POX) OWE) OX 

AFUNCF(OKVEC ECE CLERPSLNE OX) DEL TAF DOK DP OGAMMAS ) HX PD OALPHAF) UX PORE 

FAC & (HPLANCK RNAV /AHO] HODES HH ERER SIM MHES 

FAC 8b e LS TARG/SSYRTE( 7M) 

YPROR2] 

PUSTsP 

FV} 2%69 

CALL SFCONN (REGEN) 


TURNING POINTS FOUND FOR Vo g VE TASVEINGHIENS 
INTERPALATE COW) ANN OV TFTA 


jJJ2) 

J<s 

[Js 

1K 2s) 
HME See] 
“NE Ge] 

PP CTE VET) Pole? 
Pyie2e7(1) 
WRITE LAeEeTD GT 
an QQ f=le* 

LE CP OL Pere) Vezete? 
Pyztuel 
TEVOCTIJp|evi tb) 

CP percly 
wo Tu, q 
JJz IS] 

ry ( IJp2evt) 
ws(JJe2!? 

TE CPV ITD O19) AeTed& 
Tez] lo) 
arth yey 
SVEN CP ysly ( 
rey oT a 
Je zlvet 
EM Ce heey 
tel Jere? 


COMTI WE 
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9 ONWPF NE DV 


14 


a 


ee ed 
WN 


14 


15 
16 
17 


18 


19 


TF (JJ) 1795299 

CALL NTIPSR CTEM gtlQUEV oP] ttOF Se TJeJIoVNIN SVN OINSC oMREGONUST ) 
Loces 

TF CINSC) 979017997 

FO 1) TeleJJ 

ty2I4T) 

7CTIVECTT) 

CONTINUE 

IF (JK) T4014 et? 

CALL NTIPER CANMTH MAY EV EMT SHINES g fH ga JH ITM GUM Gg PNET gMRES ICT) 
LOC¢=9 

TF CINSC) O79%4997 

AN 16 Tels 

Ke IR) 

OV(TKYEAT OT) 

CONTINUE 

TF CTQHK) 1791%17 

PV LE<GFINICF(O—&) 

orrve|en 

nO 1% Tele" 

PCT IRGCEIMICE (VIP) Y+FVI 

CONTI UF J 
CQO TA 99 

INTGSINTS 

“UNE As] 

2 J 

VEVI1 =% 

UNFESel,e 

eval Vy] 

TF CEVI) 27027%97% 

CALL STORET (CV o7 GlIEV KV HOE Nh MVM aN Ge TE Ce IAT GNIS T) 
LCCcs! 

TF CIMSC) 97071 9°7 

EvV1T=FV()) 

PVT NOW COMTAING THF ZERE eAyNT Enrany ny FYTOPEAL ATION, 
Pyye7OTHN 

TF (PPM) 29995994 

Fv(y)2eOr 

SO 28 ete 

AVITP SP CIINCF EV (OT) ) 

ronTp ene 

25 TA pe 

TF (CTAMOP) PHOT rEIZRE 

Ev VFSROetiE 

an. TA 99 

TPLL MTORER CY pAV LFV EV HOE oa MeV IM ON oO EMEC POE GUT) 
EV (1) MOY CONTAING (RF. RY FXTRACOLATION 

VOCE? 

T= (CINGC) 07,78 5927 

P( VJs "6 

2017276 

eEACpPper ey) 

vite" )ev(1) 

FoUey VER Ye) 


2ITFBJIddI BAVIATBIARI,IABVIAIasIAVAgggas3z,IQIIasavvIe3sIa BBW. 


tTIIVNIIITAIA 


1297 23 


J 


>; 2 2 


Ra 


29 


71 


ap 


1% 
14 
18 


Ve, 


a9 


ag 


OVC TeDPAVel) 

Cae Te ye 

eer gy? 

TF CICBEV) 17 5%9524 
Pam Ts pe 

eVvisFry; 

ar a2 9 
At ry2%t 
vVOTpev( 
"OTe se 
AVI P petty 

OWT Pe 

TE CTOPE VY 89934037 

Ey 2e%(9) 

TE (MPOAMLTTY WME YI? 

“e] 

evilyene 

CALL MTPPE™ (Pg Vg EV GIF gHAEE gMD gM OVUM IN GUM gINEC DEA QUIET) 
tN => 

Te CENES 1) 826926097 

NOVO] 

VENMy erie’ 

VENOPSHIEY 

ZUNE y shee 

A(NDy2reF 

AOCMP eX, 


vit?) 
> ee 


yor % 


aearen ¢ 


av (sn) s 
vitje 
PCL Vers 

BY CP) SPAT) 

verge cue ye 

yNeVve ye 

ANE GS Sprnee 

erg 

WRITE (KeI*7) COVED Per (feted Ne) 
WOTTE (Ket *OV EVET Pe PIT) ofa? oo) 
TE CDE YB gh Aug BOP 

es (YE Te eC TS YS INEN GT, 

Gy OI LE oe Bog 

ATU TL SVC TASELOSTE (Tay paunen 

SEY SOT 2 9F S(T) 

WIT peGFiryrr(orwen sy 


ayo ep eour Ty 

MeVETTEMeeTy 

Gi ee fe ia 

a ae i 

clo Coat 

FALL ST TE ED CV Ge AM HEN MPa AQVETA SVE IMS INST SMREG QMHIET) 
1Aaeed 

299c tee }7) WIM EE gps nr vie MOP MITE T STAG gate 9ne 
Mere be gt tsy 

3 1927 22 'a% POEM AVL oe Pete’) 

T= (2952) Crear? 


AAR? fet ett 


Tee C Ty RAEV ITY 
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Pe 


770s 3B9N3D3 39 DDD 


72 3 3 


yO JoIdI3aD 


6®oeosa33 


tl 
11? 
1}* 
114 
11s 
TIA 
117 
1A 
lcs 
120A 
123 
12? 
123 
126 
128 
126 
127 
12° 
1290 
139° 
17} 
13? 
142 
144 
145 
124 
144 
199 
149 
14” 
54. 
142 
141 
I44 
14e 
144 
147 
4a 
149 
15a 
16) 
169 
169 
1&4 
18€ 
TEA 
197 
159 
169 
JAR 
A) 
V4? 
142 
VEL 
ye 


UCL eOVvel) 
cOMTITNUE 


INTELRATION CONETANTE 


STEP 2,9 

2192191 

¥Z2(2) 22) 

V2(3) 2A} 

1701) 2A 

XGAUS CL el peeA 1288571781297 
XGAUS (107) te 1S 16467412993) 26 
HGAUS(1 93) ee AAIVFAFFAIIVENLY HA 
XGAUS C1 oh peed *AARIBPAVEDLTE 
XGAUS 01 o8 Pelee XGAUS( 104, 
XGAUS (C106) Ble YGAUS (102) 
KGAUSCL oT) ele ¥ZAUS(1 07) 
YGAUS(C1eB)ele-¥GAUS( 101) 
AGANIS TL ol pee MENG IG7KE91481 28 
AGAUS( 1Lo2) te llTIIIAVNG) PAAKKOT 
AGAUS (CL oA ee TSBRE3I7DI07INNLE 
AGAUS CL oO dee IM1AGIAOPAOCOIAY 
AGAUIS( 1 eS) BAGCAIIS( 2 96) 

AGAUS (2 oh) BAGAIIS( 194) 

AGAUS C167) SAGAIICG 162) 

AGAUS (1 oF ) RAGAIIS(1 92) 
1712) 26 

YGAUSEC 2 ol)eeKVAIPAEIADODAN IN 
XGAUS (2 o2 Dee lAVIVEADKATAAPE A 
XGAUS (2 eA Pee *WlAQWLCANKP4ND 
¥GAUS (206) 81 e—KGAUS( 797) 
YGANSO7 oS) BLeHXGAUS(2 02) 
¥YGAUS (206) 81 ee ¥GAUS( 261) 
BGAUS EF od Pee MMEKADPIAUAIONROR 
AGAUIS (PoP) ee IANZANIAKEILAIA 
AGAIISC 2g A) eg PVIPAMELOKTOALALE 
AGAUS( 294) ZAGAIIS( 293) 

AGAUIS (2S PEPEAIIN (9927) 
AGAUS( 2eBVENGATIC( A oT) 
12(4)#4 

XGAUE (CVV Pe eAAOAAIONGMNIATVE 
VGBNE (BoP Pweg AATV*WOSPOIAIRID 
WOAHIEC Ayg Axl] ge—VAAIIT (49%) 
YGAUS (2G) By er YEAUIS (M9) 
AGAUS (Aol) tel TIWTADPEROTIIDT 
AGAUS (CF e7) ze3IPBOTAN7II“AVAF2 
AGAUIS (293) BAGAUS(3 92) 

AGANIS (6406) fAGAI'S( 297) 
7601)27603) 

WRITE (&0117) 

AN AG Te} »** 

yEMITe Fy ( 1) 

AIC DENTISTS? 

VABNISD 4) 

UF eNeNP 4} 
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164 
147 
16° 
169 
170 
171 
17? 
172 
174 
178 
174 
177 
178" 
179 
189 
18} 
1A? 
183 
1PG 
TARE 
184 
1@7 
18 
1A9 
yon 
1°1 
192 
193 
194 
yoQs6 
194 
197 
yoR 
199 
pan 
any 
on? 
Dara 
2AL 
an 
Ne 
Dn 
?aa 
ana 
1A 
21? 
212 
212 
214 
218 
214 
217 
2 yA 
210° 
220 


44 


4s 
Ak 


VAN 


4? 


48 


49 


ga 


8} 


5.2 


63 


bans 


&'5 


£6 


a7 


AN 44 Jee gue 
TE (TEMOT AVC IP) 48 ehh ohh 
CAT TNE 
T6GaNfe} 

CO TS 46 
1SGeJ 
“NOTA T SRA ICD 
VMTNE% 
FFA %e 

tet Sey 

Ltel 

176%92% 

iciplek samy 

END2%, 

1sG2% 

yeaer 


COMMENCRE FINSTNG TURNING POINTS BY INTEGRATION 


CALL SECOND (2EGTA) 
ASEVYUTNS (TEV PT AUMIM) HS TER 

TF (LT=72) 49% 94%,87 

moa"7(LT) 

vise (919/27 

TAT NST 

roene 

ASC TST) SELVATE (MH); 

oN 

TE CTOMH) POM eh TAM 

CALL “TEPSO CV soll Vel eB eM Me MIN OTS OTS OPER OUST) 
LoCc=? 

TF CINS7) 79k CO97 

TE (GNA) ©7081 651 

he 20) 

YF (7Ruwy baad Lee pane, 

CALL TOPED (MV gQVoltIEY eM Lor ariP atv o TYTN oTH oO INS pMPEG MISTY 
LOCe4 

te qrecry 7h 7A9% 

AA CE Jed g''9 

SMENVEUCT IEF VO) 

TF (REVAL ANGIE L 

NEVO SD tLe SEARTE(SFNO) 

ev ( Jypermer cyano Js) 

FAT PMI 

SC yMSOU( LPH ttEV(DPsHIEV (MA) 

FEL MSOUCD Pee AF (PP SEV (MD) 

AQ 64 JeDeve® 

SOME CIM EL HEU L IR IVEI FIT ( OH JRT) 
ASIIMERSII Sh REMI IF SHIP HEV (PH I=} ) 
aut Teka 

eer Pear eum sa, 

GS ae SSA Wwe Casey >. 

TE (L THT YP FRM H QO 

WZO%e(LT|e1)77 

re (T7973 em 07,89 
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an” 


3 


a2} 
222 
29? 
226 
2276 
226 
227 
PIR 
270 
23” 
2371 
232 
231 
224 
238 
236 
237 
23 
2a 
24n 
241 
24? 
24% 
uh 
248 
244 
au? 
24a 
24 
DEN 
2®1 
25? 
753 
PRG 
286 
?&h 
287 
DBR 
280 
PAD 
2A) 
PH? 
264 
PAL 
268 
264 
P67 
Phe 
240 
77" 
271 
27? 


274% 
274 
276. 


12993 

NGAS#I701290) 

AO 60 Jel eXVCAS 

YG IVS OCNSHRVUYVINDYGAIISCTING 9S J) OV TN 
VEV(Je¥G(J) 

CONTINUE 

“STs! 

EPGHRTE!P THeXGINGAS) 

FEPSEF ASH 

TF CFESHC1)) 979A1 942 

TF (CEPSHI PY) 92992947 

TF CIOWMK) 106 94301% 

CALL NTRPSP (CV eGeXGeYGoE PSH? NGAS oVMIN GAS INSCoMST ONUST) 
Loces : 
IF CINSC) 97064097 

YF (IZGN) 68068 o78 

TE (GOA) 690644 A 

“ST=) 

TF (TRUK) 1056679108 

CALL NTPPSR (CY eAVeUEVe®I oFEPS oft pIGAS eo MING TS oINSC oMET ONUST) 
LOC#6 

IF CINSC) ©7969%,97 

FSUM20.¢ 

GSUM#O,0 

SO 772 Jel eXGAsS 

TF CTOQHK) 796% 7C 

MEND(IJVEAVIT REYES) 

TE CRENDC IY) 0999057) 
XGCIPPAGAUSLI ZNO eS) /SSONTFIDENT(I)) 
YGCIJ)sPTCIP* KGS) 

ESUMSZFE StiMeY¥e (SS) 

GSUMsGSUM4+YG() 

CONTINUE 

TF CALINGAS)@PT(1)) 746072074 
127692] 

NGVaANSeuerypm 

CERISNSV HE Cie 

CE GD elves 

60 To 7e 

FSuUM=s0 0 

RO 77 JIel eX GAS 

PENDEOVOITPHRVYG(I) 

IF (NENA) 29,09, 7% 
KGCIPEAGCAULCT7Z9NGII/SQRTF(NEND) 
FSUMEFSUMSXG(S) 

CONTINUE 

CSUMSAT (INCAS PRE SIIM 

GD TO th 

TF (CEO GPSFE Sa] eK -6) 79 792K 
ENN=1, 

TE (SECA ICE GM eK Hb) M1 eh1—9? 

A Mele 

IF CAR SEC FFA2/FESC) aM) 86594 483 
ENCE. 

TF (CABSFUGFGP2STEG1)-=09) 845 0% o MF 
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ONODzle 

LT#t_T+1 

FEGyT|erf ors 

AEGLaneo> 

TF CENDY FIP 7,BP 
FEGSE TAF Or? 

TF (499) 9% 429,90 

TF (490) 93929592 
GEGEGEC4Gees 

YYINERS 

ITF CTOOEG) 91497591 
CALL SECOMM (CAFTFA) 
XTIQECALPTEARPE CI AY 
REGIASAFTFA 

WRITE (E9116) FROSERNTIO 
TE (LT=H$29) 47 947—08 
F2FEnerac 

GF aGEG/OAC 

OMAX( TI SSQORTE( FHF SE see yer 
QMINC LT) rOMAK( TI ar 
CALL SFCONR (AFTER) 
KTIMBAFTER ORE IN 

PE GIN ZAPTER 
WRITFOUTOUTTAOF Ho LIU oTEVPY OUI TIS PMINGIT) oRYAK(T) oXTIQSRIINGAS) 
TEMOTTSTEMPY 

WRITE (KolT5) TEMPLT IIT) oPMINGT ) oPMAN(T) o¥TIDNGPTINGAS® ) oF oGF 
CONTINUE 

YO 98 Te2lo"! 

7UT) e214) 

Y( TeV 14}) 
avV(Ty=Ry(T4}) 
CONTIANHE 

WRITE (60117) 

DE TI IPA, 

OCCT SA FY 
FOC(APSEHGILEV) 
POCCAIVERHS POF 
7OC(G)FEHE POY 

FOC (SY SHH V\’ 

PVT ( A) SEHD’ EMV 

PAC CV) ear CY 

7OCCA SKM FILL 

PAC AV eKAHOYE TEL 
WRITE (KHellA) 7CCCLOMI OI T) 
WTS 21 

no TA 964, 

WRITE (49117) 

INT4=1 

Ao TQ 96 

mw7v74527 

WITTE (&_%39) 

nN TN 94 

vXK Et 

BG l=) Jeske 


VEV(CJVETHIMSELOAT (Joly ea 


f 
a 
a 
R 
=) 
a 
R 
8 
R 
8 
] 
p 
i] 
Aa 
A 
R 
a 
R 
) 
p 
a 
R 
BR 
p 
a 
a 
A 
R 
Q 
A 
A 
ry 
A 
R 
A 
iz] 
R 
Qa 
ie} 
) 
a 
Py 
a 
R 
n 
a 
Pp 
a 
a 
a 
n 
ne 
A 
A 
A 


191 


102 


116 


118 


EVIJISEGFUNCE (UFV(I)) 
CONTINUF 

GO TN 8c 

vxXaMg 

NO 193 J2te” 

UEVC Je T’IMSFLOATE (Jel) FA 
AT(JSPZPCIINCE(UEV(S)) 
CONTINUF 

cO TN 6&2 

CALL DGUNG OXGeMEND ANGE O TEMOT) 
G69 TO 64 

NO 196 JelehGAas 

ATC JP SREUNCF(UFY(J)) 
CONTINUE 

GO TO 6R 


FORMAT (20H NUSED IM BFESRKR & 1201He) 

FORMAT (32HC THE GV) VALUFS USED BELOW ARE// 6X60 2K oe THVIXSIHGINY 
Tyo elGlFiColoFlle&))) 

FORMAT (33HC THE AVIV) VALUES USFD RELOW ARF// 6X96 ( 2X s]HVEX oe 2HAVI 
INK oS e(GlFICel oFl le 7d) 

FORMAT (6HOHDFS#IPE Le 7e7He VAINZELS eo 7eSHe VNBELE eI e THe VAFGEIS/ 9X 
1o7He NUSTSI5e7He TPSCHelSebHe Nel So4He Mel®) 

FORMAT (54HO OUTPUT VeG-eLFEVFLS AT WHICH TUPNING POINTS ARE FOUND 
L/96X eA (2X elHVIXeIHGICY)/) 

FORMAT (6(F10010F1104)) 

FQPMAT (1H1) 

FORMAT (6XeAHFFG & LPF IS eT eRe GEG & ElSe Tel AX eGHRPENIITPEN APF A,749 
TH SECONDS.) 

FORMAT (1OHC FOF V = FSele AX eh G = FUMe%eILH V/CM, RTH = FlNe 7% 
112H AND RMAX £ FIN Te? 7H ANGSTROMS. THIS REQUITED FQ, 2,94 SEC HINGE 
2e/2UK912H ALSOe BV = F17e7/20Xe49H THE KLETI' ACTION I" TEGIALS F oN 
39 G ARF FQUAL TY 2F1S,9//7/) 

FORMAT (2640 MTPPDP UMSUCCESSFUL FOR AhellHe WREN U & Flb—%) 

CORMAT (2°96 TEC PFACHIFE™ “AYET } 

FORMAT (229% CFNQ TS NAT PASITIVE) 

END 
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386 
ART 
280 
4AnQ 
190 
391 
399 
1943 
104 
206 
298 
an? 
aoa 
49a 
4A”n 
41 
6N2 
4N3 
4NG 
495 
46 
407 
4nAr 
4n9 
419 
411 
41? 
414 
414 
478 
416 
417 
ale 
439 
420 
421 
47 2- 
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SUAADITINE NTRPSR (KJoeVo eV OeVOeXMoN aK oXNIMOKNAM GING oMAEGONUSEO) 
NTAPSQZ  INTEPAECLATIA™ PROGRAM (BY SUCCFSSIVFE RANGES) SINGLF PREC. 


THF SMETHON OF LACREANGF TS USED. THE INPUT POINTS “EFM NOT AF 
FYAUALLY SPACER, THERE ARF TWO MOHES CF INPUT AVATLARLEe. THF 
FIRE T OCCURS WHEN Mos Ne THF APSCISSAF ARE THEN GFNFRATEN FROM 
XMITP eo XMAW SAND XHe AT FXIT  & NOg OF POINTS FOUNR ® AND XO 

CONTA ING THE ABSCPSCste, FM WILL PE CHECKED TO INSURE THAT THE 
QIMENSTON T° "OT FYCFFDEN, MOAINMM ££ THE DIMEFRSION OF THE xXNe YO 
APRAVS, THE SECOND “ANE OCCURS WHEN M 16 ONSITIVEs THE 
ARSCISSAF ARF THEN SSCUMED TO HAVE PFEFFER GENFQATEN PRICP TO 
ENTR2Y¥— THE STARTING INDEX FCR THE OUTPUT AR®AY “IST RE 
SPECTFIFS,. THE FIRST POINT FOUND WILL BE AT xXOOMPFER),. XH MUST 
SF GYIVEP SIPCF THE OPNINATE FOR ANY POINT CLOSER THAN eO1*#XH TO 
AN INPUT POTAT IS TEKEFR AS THAT FOR THE INPUT POINTe eNUSFED= 15 
THF NUMYAFR OF POINTS 'ISED FOR EACH INTERPOLATED POINT. #NUSEN#® 
MUST PF FVEN AND LESS THEN OR FQUAL TC <%=, MUSER 2 BOTS 
SUGGESTED. ANY POINT FARTHER THAN NUSED/2 POINTS FUN EITHEP FN 
IS FOUND FRCM THE NUSFENZ2 ON FACH SITF&. NTRPSR WILL SCALF ANP 
RESCALF IN CRAFR TO AVOIL OVERFLOW. TE THE FIRST SCALING NOFS 
NOY SHICCFEDs STRPES WILL TRY UP TO 7 TIVES MOSF REFECRE RETURNING 
IN THE FREC® MONE, (iNSC = 1) 

INPUT 15 Koel XD) oVHoNeMo(XMIM) oC XMAX) OYREG ONUSED 

OUTPUT 16 (XA eYCoINS Ce 

NSC = * TF FROGRAY WAS SUCCFSSFULe IF NOTe INSC = le 


STMENSTON XEO2°%) eo VIED © YOLINN De YOCARWA)s KICIAV® YIC2NO)» 


ywuMAa ( Pao) 


ens acnn 

WAVE KAY 

YMIN=¥ eT? 

KuRhO ce cy) 

Sev 

p2FV2e% 

pote” 

" 15 ZEMs ITF Wey KMVING WAN APE TO AF SEN, TEM ONS NOT = 1 
PEROe THEM wOC Tbe Toe Tom MILL PF TIISFN AS AMECTT ACA, 
TF (2 BWPereé 

VEX INTE (APSE (CC YMA TR SMH) Yo] 

TF TS “OPE THUAN SIMENC TA OF KO, PETRY [8 EPOAR “ANE, 
TF (MOMENTO) 9499477 

ior tisk et hie 

TE (C¥PAVRNZYIM) Feet 

wepeayeray 

wittvaevarys 

TREVs] 

AA Potter ger 

YOCTPSFLOSTE (CT oMPE) Bete Ker 

rar ree 


PO F Vole 
wTCTyey ly) 
COMT Tear 
tle ol 
WOT SRE 
AUS TORTIE ED 7 > 
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B~ID ANE Div 


a 


NSaNUST4)} 56 
NF eN-NS4} 57 
NFP SNF 4) 5A 


90 8 IsNSoNFP 


MUMA(T)20 69 
a CONTINUE 61 
TF CXC(ME —KOCMREG)) Gellol? 62 
9 vous 2 i 63 


VIF eMPEGEM I=] 
NO 10 TeMPENM2F 
KaVFerarna] 
TEMOBXO(T) 
KOCT »axa(K) 
yO(K)=2TEVP 

10 CONTINUF 
IRE’ -) 

11 TE OxTOM = T019) 12014014 

12 N22N/? 
90 13 Tele? 
KeaNn+l—-f 
TEMOsXJ(T) 
TEMS#YJI(T) 
KJCT sx stk) 
VICT  FYJIK) 


XJ(K)2TFMP an 

YI(K) STEMS 
13 CONT TINUE A? 
TREX=#1 a1 


14 TF CNFOENS) 1991%018 
15 AO 18 JatPEG Mr 
f0 17 TehSeNF 
TF (KOC SI—NICTD) 16014017 
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APPENDIX II 


COMPUTER PROGRAM FRANKON 


The program FRANKON was used at AFWL to calculate Franck-Condon factors, 
r centroids, r* centroids. A listing of this program is included in this 
appendix along with an explanation of the input data and a sample data deck. 


1. INPUT DATA 


The reduced mass of the molecule or the mass of each atom is needed. These 
can be based on either C)> or 0),. A good source for reduced masses based on 
0), is Herzberg (Ref. 31). The vibrational potential of the lower electronic 
state and the vibrational potential of the upper electronic state is read into 
the program. The value of the intemuclear separation r and the potential 
v(r) may be read in any one of several different units. The potential to be 
read in can be obtained from the output of program TURNGPT or it can be taken 
from potentials which have been published. There is also the option of 
generating an analytic potential by programming the desired expressions and 
putting these in subroutine POTGEN. 


Also needed are the dissociation energies of the upper and lower electronic 
states and the energies of the various vibrational levels for which one wishes 
to calculate Franck-Condon factors. The Franck-Condon factors can be calculated 
for any desired rotational level by specifying appropriate values for JROTL, 
JROTU. The energies of the rotationless vibrational levels are then adjusted 
by adding the appropriate rotational energy given by 


J(J+1) 

h a u 

8n* ¢ at 
° 


where the constants Be, a> Ye are read in for the lower and upper states. The 


Eg = {Be + a, (vi) + y,(wti)2t 


separation of the minimums of the two states, Te, is also read in. 


AFWL-TR- 70-108 


2. DATA DECK 


As in the case of program TURNGPT, Zare has made extensive use of comment 


cards in program FRANKON and little or no explanation is needed to explain how 
to set up the data deck. 


Card Number 


1 


(for lower state) 


4 


5 
through 


5+N(1) 


ITEST 
Il 


Title 
72H1 


IIMS, ZMAS1, 


7MAS2, 14, 
2F10.0 


LIRA, LIEN, 
N( IKON) 


14 


XI, YI 
Format in columns 
13-72 of Card 4 


144 


1 in column one if problem is to 
follow 


carriage control in column 1, 
name of problem in columns 1-72 


IIMS = 1] mass based on C = 12 
TIMS = 2 mass based on 0 = 16 


ZMAS1 mass of first atom (reduced 
mass of molecule) 


ZMAS2 mass of second atom (or 


zero when using reduced 
mass) 


Se 


IIRA = units used for internuclear 
separation 


IIRA = 1 atomic units (a,) 
LIRA = 2 Angstrom units (A) 
IIEN = units used for energy 
IIEN = 1 atomic units (hartree) 


IIEN 


2 1/em (kayser) 
IIEN = 3 electron volt 


N(IKON) = number of data points 
for potential to be 
read in 


XI - internuclear separation in 
units at IIRA 


YI - energy of corresponding 
point of potential in wits 
of LIEN 
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Card Number 


6+N(1) 


74N (1) 


8+N (1) 


9+N (1) 


(for upper state) 
10+N (1) 


114N (1) 
through 
114+N (1)+N (2) 


124+N(1) 
+N (2) 


134N(1) 
+N (2) 


144N€1) 
+N (2) 


15+N (1) 
+N (2) 


16+N (1) 
+N(2) 


NDE, DE(IKON) 
14, F10.0 


XMIN, XMAX, 
XH, IIRA 
3F100, 14 


NL, IIEN, NET 
314 

Formats for next 
card 


KV(I, IKON), 
ETRIAL (I, IKON) 


IIRA, IIEN, 
N (IKON) 
same 


XI, YI 


NDE, DE(IKON) 


NL (IKON) 
IIEN1, NET 


KV (I, IKON) 
ETRIAL (I, IKON) 


NI, NS, IPSIQ, 
MAXIT, EPS 
414, E10.0 


145 


aint te 


NDE = 0 if zero point of the 
potential is at R=@ 


NDE = 1 if zero point of the 
potential is at R = R, 


XMIN, XMAX the minimum and maxi- 
mum value of internuclear separa- 
tion at which the wave functions 
were calculated, with mesh spacing 
of XH all in wits of IIRA 


NI. - number of expected energy 
levels with units of IIEN 


NET = 0 zero of potential at 
R = @ 


NET = 1 zero of potential at 
R#R 
e 


KV - the vibrational level 


ETRIAL - energy of the vibrational 
level KV(IKON) in number 


Same as for lower state 


Potential for upper state; same 
description as lower state N(2) 
in number 


Same as for lower state 
Same as for lower state 
Same as for lower state 
Same as for lower state 


NI - 1 prints iterations, other- 
wise not 
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Card Number 


Card 16 (cont'd) 


L74N (1) JROTL, JROTU 
+N(2) 213 
18+N(1) BE1, AE1, GEl, 
+N(2) BE2, AE2, GE2 
6E10.3 
194N (1) TE 
+N (2) E10.3 
204+N(1) NPOT 
+N (2) 13 
21+N(1) Blank 
+N (2) 
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NS = 1 prints wavefunctions every 
IPSIQ points 


NS = 33 prints eigenvalues and 
node count 


MAXIT - maximum number of times 
SCHR will try and satisfy 
convergence criterion 


EPS = convergence criterion 


lower a.d upper rotational levels 
respectively 


rotational constants for lower 
and upper states 


TE - relative difference in 
energy (units of 1/cm) of 
the lower and upper states 


potential printed every NPOT 
points 


= 


ame 


Ni: 


’ 


SSVI DOT 


in Ys ee 


= ANNWDAADTDANVAIVANADAANA 


TEV CY 


+ 


bea ee i 


NOT REPRODUCIBLE 


PONVADAH FR As Peery SVT MIT M1 He psy 

FOP e AN CALCULATES WINOATLOANWAL TOAME PTE Ae MORAAMY] trtre END 
PLATT “ALEC UL eg, 1? 165 @ASFN NP B PRAAPAS! NETAPSEN FNM 9, 
FARMS. Jur DOINTe PFET NPG Tur PCT TEAL £66 SF AR fe TT Tee MIAGNH AP 
A TOL TR CVT Te EEN TR FMTEOMN ATE FRO TURE TA ye YMNTELOATING 
F.2T%¢ jue yen AOA VE JRE DOTENTTAL te Turn FryTrpProry *TrFn ov *# 
Fuycyv yer Ne kHeE FAP As/PawE (Oy, TUE AITO? oMPMeU pe CY TOFOAL ATES 
AV A FIeeTine AE Jur FAN NFAA/OBSE (DY), TNE GA NFELVEN DATE TAT 
PT. TON EH THE SOHOANIMGER ERATION nyo Te SNLVeEs Ov eoMnOnITyer 
GCHIA® THE vAVEFSINCTIANS ADE TN TEGPATEN Py VEANE AE fo cyeneAm 
BIE Teer ee Fe AVE AD TH TESCO AL, 


Syohit weed. Wr a Fab ag CR Tyr Tur OATENT? Apetres ar AN ADL! 
BEL ae NEO AE LEVELS ee Ty DE Fe AM ycde 


roles. VEC VeKT LAM AVAL PIA Der EMAMAAD SK LVAOHO) CET AL CAMINITA TS 
PPOSIEM ETDS TT NSEZD OPIMSE PMD e7IPALA PV eFIENLAD oECALT IAA) g YPM (CE, NVELAA 
PPOVOAL CHD TITER ITE LAD eV IEM TEVA) oP VELAN) SPECI TING WAY OP OMNTONI AME TAY 
Be VAs ( Fo 7 ePMAMCN( We) POL V( WAND HV IID MIP OPL IMIS IFI7D QDOF 
AUT CD VD OME TEV) PPLE) oe SIMPOE IF Pe LIAD TITLE LA) aver gQiC TA gO TS ge 
Be XMIEN VAM GDH MET MITE PT oP CANIT GEAR Gg THe TMM EME SIDE gh rhs 
AMT DION 

NTMENG TO! KE TATE EOD) 


THE FALLOUT COMMER CARRE Dro COTRE THE DOFDADAT IAS AF KATA, 


FIRST CARD TH MATA RE A OME TM TOL PT TE A PRNALEY FALLAWS. 
VERY LACT CASON Yer MAT® Sse? OF A OLANY CRs eee 

myre 2 a es Be ial] ny De PSE Se 
MP yy? CARSN TR RATE ke eer Ag AAA Ee * 
CARDIAC “AM IN| TQ re ET re dS cecnr 
MESS OF Taal TH: MTN er, 

TE PNG 7s MAET TE TeV EM ME THE PEQUCEN “AcE fe Pur Pye cwopee 
TTR 2 hs “MARE HT TS eer NACEM AN 71> H 99, 

tres ze Fe HAO FITTS ANE MAGEN CU AVE 2s Ye, 


THC ALIe 1a 2, HERE 
: 


CRON wae VPaee vasn 


CTO OTS FO SRT Be QEy Pie AIGe gs S OEE CONICS Te BN A TOUT AN IRAE ES 
DEAN Vag SoH SY 


{F CTTEST) Ferree? 


TECTR a CMEC SOON E LER eens A RST eG 2 sure me ryan | FT CALC LSTe 
ra 89 SBFanR rent enor 3 Tov TM ]er PA OAL Aree Ore eRe TAR pte 
waar 


OT E2 2A. WET LE 
BO ATY Whe PLOT. 

CMTE CT) = 2 WIE Paven. cA Fr. CF ATE (SO) om Po Tur Annee ey 4F¢ 
HEA PI, VERN TEI 

BFAD TR, vOTe re 7, 

DEAN Pre TTP MAG PAD 


COTA IN FAD NPAT TY. PATENT HPL 


FUN Ae MPA cero 


“ury? rere nm apr a rn a ei YES SC VANAV a VES S TOE ® 7) ated a Cn ia melee S = 
os eres ‘ ‘ 


>» ep» re oe 


Dp 


ad 


DPp7~ DPepre PPD yY> > 


‘> 7 » 


~ > 


> 


D 


Ss 


- 
73s FQ 7 AIPAAN 


me ee 


LVN 


bY 


AO NAAR EV VANE SY OVA Oe. 


Ds a le Paces ON <a 


“E™S CAND Ter MTA ANA MAM ATE TANCES ANT THE SDACTON wet, 
MEAN AF, VETTE AV TE VHT eg] IDAT 
SEY PT Tere CAPNS ADE THR VALTIES CE J ENQ Tue LOVED ARN [INREO STATE 
MOAN AT, IMAM § JONTH 
NEYT CAON ThE PATE WAS TITRA. LTTE Ns ANA THE StIMREFO ANE ONTHTS, FACH 


TN 16 FORMAT, AL VELL AS THE FORMAT STATEMENT waetcCH CASTROLS THE 
REASTNA SE TEE OATITS CIM COLUMNS LIRATPPemEHD FKAMOLFe (4FIHeA) 


CTF 2 Ye FMEOMV TH J ATHMTIC LINTTS, CY Belts = M%eIMH) FoeVeo) 
TY 2 %, FMNEORVY Te Th LT sCeey 

ITE fF, CMEOGY Te Ty FLECTONM VOLTS. € 1 Fee © "VAR AM T/C) 
TTON 2 Ve TTCTAMCE TS TN ATOMIC tINETSS (€ 2 Aalto © eS 291HK ANTS) 
TI2A & Pe AISTANCE TS TN AMAETR ANNE, 


FOLLOWING CARAS COMNTAIM THE INPUT POTENTIAL POTTS. 

IF YT 8% MESPTPEN T9 ISE A EFLINCTICON TO CEMIFDATE THE wunLe 
POTENTTAL CLOVE, THE EIRST 12 COLUMNS CF THF TIPAg FTC. CAON 
mie T OF MLAS. A HOLLER ETH TEXT MIIST FE MrIMCHER FP 17S FD. AS 
TT TLE TF OP PTEN, THE SNF XT CAPR TN Sree A CAGE CANTAING 
WTO VAV OTC 


HPRPED STATE YTS PFAN IN FIP STs LOWER STATE LAST 
THE SAME VALIIE OF TIEN STIR ANNE e XHeKMMIMOXMAN OLIRAY gNET 
el TEN] 8G y AE MICEN FAD AAT THE GOOHINA AMM HINMED CTATE 
wwAN 29 
OFAN Arq TIMP STIEM MI ITHAND oS (MN TEQ*T(T) ofl eI) 
ver Tere weep ee reese o 
vAVeY AY TES 
Via yriTe 
apenas ( Pumery 
MEAN NTEOMT Sey PT (Tp eVP(TpeleleIV ON) 
mNAer( Tenens 
APDe Arad 
MEAN DTEOMT.L VI (MIDE SYTINIDY BNNVEN®N? 
f CUS Bra Gh SU dd oo 
BAO We) ett) 
Tee yaet el at 
gees fore? 
BEAQO NFER YET. YEE TD pevV CT r peavV El Prev Pl Tey 
VEL av epee 4 


PAT pee 


NEY TRON CRT ATE NIE AMR TUE RTE CACTAT IAN EMCO nY wERr ane fe 
JOIN YE PEON AT AE MATFUTTAL CIPVE T& AT PeINFIE TY, AD AWE TE AT 
PaPEWTE POMP es, ATCCACTATIASN CREOTY WAT CAME rE PE AS AATENTTAL » 


OFAN O21, SOE ENO CT KAA) 


PVT CROP YAS 8. ANN PPE CE FYPECTED FNEPGY LFVEL Se Apo MET 
SMTTH LA TUE SCARE FNTFROOETATIO*N FOP FYPTALIT) AS MAE SNTS FAP 
vue NATENT Tay CHIOV APNE, 

Tv ALTO COM TPT Tee VAOTAMLE FAPMAT END GCFANTEGA FNEORV LFveL © 


148 


7p >> > DP PPP PFPePeVWe FR PPr PPRPTPrPYPPYePrPPrrrZ PPP P DP ?DWP? 77> DP 


DpPRPn» »> PP 2P> > 7? 2 > DP 


~~ D> 


54 
&? 
Sa 
gq 
6" 
4) 
4? 
aa 
KG 
6& 
6A 
4&7? 
4a 
49 
7A 
1 
7? 
71 
an 
74 
TA 
77 
7a 
79 
An 
| 
ad 
RA 
m4 
Ae 
an 
ay 
na 
ana 
an 
97} 
.? 
oa 
96 
ge 
96 
9 
na 
an 
yan 
1} 
19? 
qna 
104 
pre 
14, 
wn 
ina 
yAN 
TA 


aN 


NVAANDAAAIADVAAWDAGD 


NANDWY 


AVNWDANY 


ANN N 


OEAN Phe MRED OM SDPO eMIET O(N EMTET PD efelele) 
HE Penery Cpe rer 


MEAN TM EMNFORVM LEVELS 
PFAN NENT, CV VIT POF TOTAL (T Dela? MILIAN) 


REAN 7% Teer CYDEOMMENTAL nv 0 TALLEY Ove 
PEAN AR, (AVEC TV eTeLeMLINNM) 


AAING FCYPST FAD THE TIMMER FC TATE AMA THEA THE LOWER CTA TE SeCUDAN TANG 
FO COMIAT IN IGT THE FALE Re ME THAR 

AD 8h SCHP MARTE TEN END TUF PAASIANASIOEN:, Ye NeEeeprnen = ye 

Me Ne 7APF AMD te He CASHIN, HCPL me TAAAY_L 19412 

EOP THYS PRAGOAS® THE CANE TANTO ADE 11CEN Ae FOLLG IC 

TE NT = Ve PPINTS TTEPATIONS, ATHERWPEE HAT a 

TF SIC = Pe OR IMTS WAVEEIIMNCT IAM EVERY TOTEM DOINTE, ATHEDWEEE NOT 
TE MS & 3% PRINTS FIREMNVAL TIES AMM MANE CHIT 

EPS 1S THE CONVERGENCE CPLTER TOs 

FPS TS UM Ter CAME TET Ae TE MATeEMT TAL FOE 

MAXIT 1S THE SAY] PIM MIRED AF TEMES CCHD ITLL TRY TO CATITSCEY TRE 
COMYEOREMCE COPTTEN PAS 


LEP rne 
VLIMa? 
PFAN Phe MTN eS TOSTQ eS MAVIT EPS 


MEYT CAQDN COMTALTMNS MONT, 
THE POTENT TAL 1S DOTNTEN AY FVERY MDAT HALT, 
MPAT & ER TS CUICGFETEN TH Tur vicerD, 


PEAN AY, NAAT 


THIS TEOIPATES COMMENT CAP AE ON Tur PIFOROATION AE NATE, 
TE NA FHP TUED OPARLENS COLE Ag PEMFMRED TA ANN A ALAM FAO TA TE 
RATA SEP OY 


ANTOA( TV SAUA SIlg 

ATPA( PP RAUAUIOET, 

ATESN( VY EKIUIP gt'g 

ATEN Ay eR Oe RS 
ATEN ( FZ) RAUF OV 

ATS) peer) F299 

AIMS C7 SANOVAS TA 


AMUITO NT FE Pee tev AIMEE ALIN AM Re rNAee. 
PRTWT VP ANT IC 


PRINT Vo, TITLE 


DOTNT AY_ JRATY gyrate 
POPHT The wae ere pen TrID verte 


TE (7PMATTD Rahek 


149 


> PPP r ry e>DRPr > »> DP > PDP D> 


p> 


2 Pe ?F PrP TBP VWPDO KRPVFRZVPTZ PDRWATVZDZ DF 


>_> => > Pl D> 


v1 
TD 
1 li) 
114 
176 
VIA 
217 
17¢ 
119 
124 
171 
17? 
194 
24 
y76 
194 
he 
178 
12a 
yAn 
141 
949 
142 
i WA 
Ye 
134 
127 
yaa 
149 
yan 
147 
14? 
142 
You 
146 
Yue 
147 
%44 
wan 
TRA 
vat 
169 
364 
Yeu 
eR 
.h4 
wa 
180 
769 
yan 
%#) 
149 
Vat 
TANK 
ee 


NVANWDPY 


io ee ee) 


i 


CQINT Ae ATPCUCTIMS) o74AST op 2MASP 
Go tT. 7 

PRINT 8%, DEMSCTING 2 e748) 
CONTINUE 


PRINT THE INPUT POTE" TIAL ® THE DISSOCIATION FNFRGV>e AND UNITS. 
TE PEAN IM, 


PRINT 80%— ATOATTIPAD SVTENMIT TFN) 
MEY ENC THA) 

PRINT OMe (YTCTIOYTCT ef @leNIMON) 
PRINT Ae NEWT EAN 

PRINT OF 6 ViNTN e¥MAK KH OMIOA(TIRAT) 


PRINT THE TPTAL FNERPAY LEVELS ANID THEIR IMI TS, 
PRINT 9% NTENCTIEMNT PO (KVOT DoF TRIAL CT) of 22 oNLIKON) 
PRINT THE CONVERGENCE CRITERION . 
PRINT Ag FOS SNTFMC TIE MN) 


POTAT WANE 
PRINT Fo, TITLE 
PLOTTING GENFEATEN POTENTIALS coe 


TF (7MASP) M—eAeO 

P°*Upe PAC? 

An TA 14 
PMtPePMACT@IMACI/( PMASTOIMATIY 
PIMC DV) =1,” 
PIME( PE SII*ATOL 

PIRACY VE1—% 
PIOA( AP) HY eRONTAL 

Perse 7OUIM7Z TMCCT IMG) 

TIFNC DDS 2 AKU IBAOF? 
PIENC PY e7MIFTLAINADEE RD 
PIEME BAP HILL TAINDVIEED 
FCACMBAN SG ONVIINDS IM 
M>eN(2) 

ML AeNL ET Y 

INMeNy (Dery 

AV 1} Teteltin 
VECTIEYTCTPOZTEMCT TEA) 
YICTIVEVICTISZIPAITIOAY 
COM PME 

RESNEOZTEN( 7 Perey 

NEC PVENFCPVE7ZTENLT TEM) 
WY TNS KMINGZEPACTIRAD) 
YVAN SYMAYHZIPAITIRAY) 
Whiz KHt® 7TPA(TIRAT) 
EPSSEPSH7TEM(TIEM) 

TNAL SLOT) 

HAY? Fey ofst0v 

EYOTAL (TSF TPTALICTIS7IOE MCT TENT) 


>pPrPP SYP D> PP PPPP> PP? PP PPPPePPP?> YP PRPPPPPP PPP PP PPrPY> PrP PP? >> PPP PrP PrP YP pp rp 


164 
167 
168 
140 
170 
171 
179 
17973 


174° 


178 
174 
177 
178 
179 
189 
1A} 
182 
18% 
184 
148 
184A 
Vay? 
1ARe 
yaa 
199 
19) 
192 
1943 
194 
yo8 
194 
197 
1OR 
yao 
JAA 
ni 
an? 
ana 
ING 
ans 
anK 
In? 
9”Nna 
ann 
21”9 
11 
21? 
9113 
214 
218 
214 
217 
210 
919 
22a 


1? CONT INE 
CALL PATENT 

c CHEC™ FOR RESCALING NE POTENTIAL TN ZERO AT Re INFINITY 
TF (MME) 918512 

13 AO 16 THlem 
VOT EVETIMENEC TEAM) 

14 CONTIMUIE 

15 TF (MET) VretMelh& 

16 "0 17 TeleYMAL 
ETRIALCTISFETPTALITIONECTKON) 

17 CONTIME 

18 NO 19 Tele TAAL 


FCALC(T Ie TRIAL(T) 233 

19 CONTIMIE 246 
€ 228 
Cc PRINT THE POTENTIAL CURVE GENFRATFD (ANGSTROMS AND 1/0) 236 
Cc 272 MRITE CL eI STATE oT KONEM 247 
c RUFFER NUTECT el CXOCLD oXO(M)) 23a 
¢ 138 CONTINUE 239 
c TE CUNT T 0191980 76097400 24N aan 
Cc 290 RUFFER MITCLoTIOVEL) eV OM) ) Pu} 
PRINT 98 24? 
NPONTESNPNT #4 24% 

a NPOTS BNPOT HS 246 


NO 26 Tal eMeMPOTS 
TF (M=TONPONTG) 297107) 
20 JFINS (MoT) /NPOT 
TF CJFIMD 279979299 
a1 JFINGS 
2? AO 2% Jel eJFIN 
TORNSTENOQTH( Jal) 
KPRN(JISXN (TORN HB INTAKE 
VPRN( J) a(VITPPND ENE (IKON) POFACY 
23 CONTINUE 
TF (VPAN(TI—1LeFA) 28 e2h eh 
24 PRINT Se (CXPRN( J) eVEPM( Ste Ja) oJFIMN) 9 PORN 
G40 TN 24 
25 PRINT 275 C¥PRA( JS) pVOPMN( J) eo Jat oJFIM) » FORM 
26 CONTI ANIE 


a7 MTEMQ eM 26 
LLken 261 
VLKeN 26? 

c FINN THE FNFRGY LFVFLS THROUGH USE AF SCUR, 265 


NO 4A Tal oMILIVON 
CALL SCHPA (KVIT) SFCALC(T)) 
TF (SCHQ16) 2%97R071 
PALL ALEK 4} 
TE (CLUK@LL IM) 29570,73 


29 MeV TEMP) 249 

NO AN JaleM 27A 

. an YOCS) aC Ja1LPexXHexM IN 7! 
AVSUMESCL PEASY N( 182 GL OS (2 VEHISYA( 2 )HEDSE (MI BED/YN (MY OED 27? 

VSTMD ara] 274 


AD 27 Je2eKSTMN GD? 
AVSIM re HE(IDEHI/XO( S)HO24H 0S (S41) FH2/V0( Jo} ) PEDEARYVSHIM 


PPrPPprPPrpePpPrPrPPpPprpPpPp>»>PpPppeprprrrpprrrrpeperPrrrrrprrrprpyrprprppprPpPprpr»pp p bp ppp 
N 
e 
B 


bl 


ae) NANANANs 
~ 


N 


» 
~ 


CONTINUE 

ALSO WRITE WAVEFUNCTION ONTO TAPF 
WRTTECLITUOMNS MTOM ST oh! 
TKOM<POTFENTIAL eMILIKAN@MIIMAFR OF VIRRATIONAL LEVELS FOR 
POTENTIAL IKOMeT<NUMREP OF THE LEVEL eMeNUMRER OF POTNTS AT WHICH 
WAVEFUNCTION ITS TARLIILATEN 

AYEFER MIT CLeDIECSOCT DoS EMD) 

ANY TA (3%o779— THON 

OUT WAVEFUNCTION INTO EXTENDED CORE FOR LATER USE 
INFCwlo(l—1) ar 

CALL WRTTEC (SOLD Oe IMECT OM) 

a0 TA 30 

NO 7 TheleML? 

INE Cae lome( remy) 

CALL REARFC (SIHNCDTV oT MET MY) 

NY AB Jato"! 
KOCUPSCFLOATE( Jel OKHoUMEMp eS) 
CONTIANTIE 

CALL SIMP (SeStieMeKHeRESULT UT oIK)) 
CALL SIMP (SexXDeMeXHePCMNTPA (ET oIK)) 
AeXMIN 

AO AR Jat at 

VOCS BAPKO(S) 

AsAeyH 

CONTINUE 

CALL SIMP (SeoXMeMeXHeoRPCN(T of KI) 
AgXxMIN 

NA 2H Jay oh 

KOC IP SABKOIS) 

AwzAeXKH 

CONTINHE 

CALL SIP (SexXNe eM ePORCMH (TSI) ) 
Aaxthy 

NO 27 Jzlo" 

KOC IP SABXO(N) 

ARASH 

CONTIANIIE 

CALL STYP (SexXM Me KH eQOPRN(T AIK) ) 
CONTINUE 

AVE (T)sFACMORY SIMaXHs 4, 
FCALCULIENE(T KOMP SFACHUSECALC(T) 
PRINT 9%, FCALCITPeAVETT) 

CONTI IE 

OQ, vo, TITLE 

ONTNF AA, LLY 

OPIMNT 1% YSFATE (TVA) 

OP INT 47 

AY 47. TeremL Tan 
FTRIALCTIS(FIRIALCT PONE CI KOM) POF ACH 
NTFFSFCALC(T)<FTRPLAL(T) 

TF CT H9) 6904206) 
AGTSETPTALIT IF TRIAL (f—1) 
NGCBEECALCITYOFCALCI T9212) 
ATENE SRA CANLT 

DOINT JTF, ALT ANAC ANIENRE 


PPP PRP rPPrPPrPPPPrPePr?PrPrPPPpPYrPrPrrrPrPrperrPrpr»rrrpPrrYrrprrprpPerPrrprrrprprr,> PrP PrPyrYrprprpypyrrprpp 


42 
re) 


nh 
45 


46 


47 


4A 


49 


§1 


52 


63 


ORINT 1920 FVOTPCFTRIAL CTI oECALCIT) eNIEF 
CONT TPE 
PRINT 1%ho YSTATFE IC TENN) 
PRINT 106 
AN Gh Teoh reo 
ORPINT 1NHe KVITPOPVE ITP AMVE IT) 
CONTINUE 
F0 TN (HAF) TYAN 
KON! 
an TO 4 
CONTINUE 
AT) 47 Jel ofl 
AO 47 Tele ? 
OF SRP SILT (Se TY) 
Nah 6 I9164K 
BCNTOEN( Se IT hVEPCNTPN( Je TV PBA SOFS 
MAREN LE IN1KA 
PREN( Se TMP=ePOCH( Je TK VEANSPES 
RBREN GRIN AA 
PRRCN( Joh EPPPCHN( USI KYSA/SQFS 
ASHEN eKINIKA 
PRPOY( Jel” pePONrn JePK YHA /OES 
CPESUL TI So TK) 2OFCHPFS 
CONTIMNE 
PRINT 197, TITLE 
NLISSISMLIF LIN 
CONTINUE 
PRINT 19%, CHV JD) oJeML IS oMLIE) 
AN 49 Jato? 
verJjq-} 
PRINT TMs Kel(OFSULTOTL& oJ oTLO ELIS eNtLIE) 
TE CMON( Se VeNF EMI FN TN 49 
PRINT 124 
CONTINIE 
PRINT 176 
TF (NLIF OOF eft) GN TH £4 
MLIS@NLIFSISML ITF eat 15940 
TF (NLIF SL TOMLI FO TA 4f 
NE PE RMLEGATAI? 
CONTINUE 
Sump FCF 
MT ST Jel oftt 
AC &Y Set e'L? 
SUIMOD (Dee SEMPP (C yPsPE CHIL TC Ie) 
AN §F Jel eML? 
AN 89 Jaret 
SIIMP (CIID SSID JIPSOPCIIL TE Ped) 
PRINT: 7a, TITLE 
AO 63 JeleML 
IMI sJq-)} 
PRINT 11% JM) eSUMPMO( 3) 
SUMPP (PEN, 
AD 8A Sel elL? 
SIM Le IIa? 
PRINT LPLe SIMD eSUMOC IS) 


153 


7 >> pPPpnpPrprpe» PRPPP PRP PPPrYrPpPPPrPPPrPPF PP > PP PPPPPPPP> PY P> PPP > P DP RP > > PDP 


avy 
a29 
133 
a. 
are, 
ade, 
aay 
ata 
aan 
Aan 
aay 
34? 
442 
ay) 
An 
246, 
way 
24R 
2409 
259 
a6) 
aaa 
464 
254 
age 
156, 
157 
269 
aga 
aan 
161 
aA? 
Ahr 
he 
cy a4 
764, 
AA? 
ane 
aan 
27a 
avy 
a79 
a72 
atu 
ave 
a%6 
a77 
a7a 
aro 
aga 
aay 
aang 
aaa 
aay 
aae 


| 
| 
| 
| 


4 CMO Vel, A @Ag 
OPHI'T LIA?) TITLE A 287 
ML VSeVe tt PEslA A 32Aa 
ag CONT TA a ana 
PRINT 197%e (KVE JSP eJSeMLISeMNLIF) A aQn 
AY SA JeloML? a any 
2 Jo} A 209 
PRINT 1146 Keo(RONTROCILM et) oILKBNLISeNLIF) A 393 
TF CMOM( eS VeM& eo) GN TO 86 Aa 39%. 
PRINT Y04 A 208 
&6 CONTIMIUE A 394 
PRINT VIE A 1209 
TF CNLIF SCENE Y CN TA &® A ane 
NLISRNLIF 2 CPL Teeth 14472 4 490 
TE CNLIF OL ToL GO TH BH A GAA 
ME VES EAATARNSY A 4} 
4 CONTINUE A an? 
PRINT LlH> TITLE A 403 
MLISSISMNLIF SVS A 44 
SR CONTINUE A ane 
PRINT 177% (KVOISD oJeMLISoOMLIF) A ans 
AN 59 Jay oft? A 46n? 
va j-} A ane 
PRINT 1146 Ke(PRCNCTLW oD) oILKENLISeNLIF) A 4no " 
TF (*ONG Dee eT} GA TO 89 A 41N 
PRINT 116 A 41) 
&Q CONTINUE A 41? 
OPINT 1156 A 613 
TE (MLIE CGF ONLY GO TO AN A 614 
NLISSNLIFOTSEML TE HAL 15494 A als u 
TF CMLIF OL Teh) GN TA £A A al& 
MLTF SMLeAntTrannr? A 41? 
aN COMT TIE A 438 
OMIT YI%.e TITLE A 4619 
MLYSs PE TF 2] A 42” 
41 COMT Pe A 4?) 
ORIMNT YI7%e CHV I) oJeNLISONLIF) A 429 
AA HP Jad ott ? A 494 
vsJ-} A 494 
PRINT 112%. Kel PPPCNI TILK oJ) oILUSILISeILIF) A 428 
TE CMANG eB Ve H OM BAN TAN 49 A 496 
PRINT V4 A 427 
4&2? CONT TSE A 4?a 
OR TNT 116 A 429 
TE CMLIT SCO Y AN TA A A 432A 
NLISS!L ITE SVEny Tee IL 1S477 A 4%} 
TE (MLE oe Tet? GA TA 4 A 41? 
PE VE SNL ELATARRY A 432 
43 CONT TANUE A 416 
PQW7 11e@, TITLE A 438 
Ct] i ee eS Oe A &%% 
44 COMT PME A 427 
PRINT T27%e CKVEI JID OeISNLISOMLIF) A 4&3" 
AN AG Jad ght? A 429 
va fo) A 44An 


45 


a) 


67 


AR 


49 


PRINT T14¢ Kel(PPOPN( TEKS IP SILMEML IS SMLIE) 
TF (MONT DAS eNE CNY FN TN &5 

POINT 194 

CONTIMHE 

PRINT 1146 

TF CNLIF OCF eMLI GN TA 44 
MLYISANLIF ST S"L PF el 15473 

TF (NLIF ol TeNLY GO TX KA 

NLIFENLSGNTOSS§& 

CONTINUE 

CUNCH 199, TITLE 

OUNCH YIL—e MLO ML 

DIINCH 79° 

PUNCH LIF, COT sJePF SILT Sol po t=) oh peletont 7) 
oOryecy 1237 

PUNCH 1249 (UL eJeRCNTRO( Sold oteloML) eoleletit?) 
PUNCH 176 

PUNCH 1266 (CPeJeRPCN( Jeol) oJeloNL) oTatoNt?) 
BUNCH 1276 

PUNCH 12ho (OT etJeRRRACM( Jol) eJzloNL) oT al eNl 2) 
NO A? Ther eML? 

NO 467 JeloMl. 

PESULT (So TKYBOCNTON( JoTu)/PPCN(I,IK) 
PCNTPD( Je TX) B2OCH (Ie Th) /SRRPON( Jy TK) 

RREN( DoT KP ePPPCH( Je TK) /PROPN( JeTK) 
CONTINUE 

ORINT 19M» TITLE 

NLISSISMNLIF 214 

CONT IME 

PRINT L27%e CHV S) oteML IC oMLIF) 

AO 49 JeteMt 2 

Wa Jae) 

PRINT 1295 (PE SULT SoTL) eo ILM EL 1SeL IF) 
PRINT 13%e Kol PCNTRA (Je I LHD oILMEMLISONLIF } 
PRINT 19496 CROCH( DoTLK Do ILM ONL IS oS NIF) 

TF (MONG Deh) eM et) AM TN 49 

PRINT 114 

CONTINUE 

OR@INT 1756 

TE CNLIF SEE OMY GA TN FA 
MLISSNLIFSISML VEEL 15499 

TF CNLIFeL Te LY CM TN AA 

MLIF RILECKIANANY 

CONTINUE 

AA Th } 

viest Le} 

TE (VL Mert Ty) 99,979,979 

DPINT.]A%—_ rE 

an Tt. } 

PRINT LAM CLI Ve! 

an TN 1 

ORINT Be "AE 

RO TS 7 

FALL EVI'T 


A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
4 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 


76 
7 
7a 
79 
ar 
a} 
a? 
aa 
Q4 
a5 
a6 
a? 


ap 
ay 
an 


91 
9? 


q3 


94 
9§ 


96 
a7 


oa 


ag 
rae 


1 


VA? 
yea 
1% 
1G8 
Liars 
yar 
ya 
jag 
nhl 
V1 
1h? 
113 
4 
VTA 
114 
117 


EFQR“AT (TT) 

ENRMAT (AY) 

EQQ*AT (ORV 

FOOMAT (TH ePFI%) 

FNP“MAT (AED ATG) 

ENOMAT (D7 2y 

ENRMAT (ATOe IAB) 

ENOURT (Teel INAS) 

FAQMAT (ATA eICAK) 

FNOQMAT (MEINAA) 

FOOMAT (OTM 0F1%—%) 

FORMAT (///4%4 THE “MACSEE OF THE Trt ATONE, RACED AM AK ebHe ARF F 
VIVA eBH ANN FYOLAS/S/) 

FORMAT (/7/446H THE PFNUICFD YASS OF THE THQ ATOMS, RASFD ON AfeSH, 
Tis ©1%4///) 

FORMAT (35H THE INPIIT POTENTIAL POINTSs R TN AfeolSH AND FNFERGY IN 
1 ABs lBHARF GIVEN PELOM, S/IUK SL AXeTHROLIX oLHVeGX)I//) 

FORMAT CLIVE I% eB oF le ebeK IM ebok LIeboF lNeboFleae4oFlNeboeFlPeboFl Nebo?! 
1904) 

FORMAT (&47HS NISSOCTATION FNFRGY IN SAME UNITS AS ABOVE ISE2NR) 

FORMAT (7/7/94 PMIN = Fle 7e%Hs RYAX & FINE Te12He SPACING & FIN27 
LeQ4. ALL TM AB) 

FORMAT (////29% THE TOTAL FNFRGY LEVELS IN A6el@HARF GIVEN RELOW. 
T// 1X oS (SHLEVEL 0 AX eGHEMNERGY OAK I/SSI1H TOolPF IAT elboFlbeIlGeFl be Tele 
PoF Ae toelhoF lhe 7d) 

FORMAT (////643H CONVFEPGFNCE CRITERION 15 FRROR LFSS THAM FQe27 e2KA 
14) ‘ 

FORMAT (374% THE POTENTIAL FUNCTION GENFRATED IS///1K oS (AM eGHRI ADs 
TAX eS PHV ET SOM) PY)Y 

FOOMAT (CIXIPF I GHelOF IAG eh (MPF IL Coho IPFI%g&) oaXIS) 

FORMAT CLXFI% Bek Lr beF IM eheFk ll rohoF IN eh oF] PeboFlNehoFl PeboFlNehoK} 
Peek VIS) 

FORMAT (39HD SCHR FINNS FNFRAGRY LEVEL G © Fi11l)eSe1SH 1/C% ANN AV w F 
TVT)e7%eS&H T/T) : 

FORMAT (29% PONGPA” SHICCESSFille (MAXTT REACHFN 1997H TIMES)) 

FORMAT (72H) THE RELOW 16 FOR THE AlNe6H STATE) 

FOPMAT (S4HOVIRS NOe Givens FNE NAV CALT. ENEORY NTFFERFNCEY 
LAX eA®HGTIVE MFLTA G CALC. NELTA G RTSEFRENCE /G4X% o 29KHOTEFFRENC 
Pe 2 CALSCe — FYVENSS) 

FOOYAT (AUVIFTE,E) 

FORMAT (19,7F VERY 

FORMAT (79K) Tur AELAY TS FOR THE STATF AYN) 

FORMAT (VEAP AM AUT IS Mee AX QQHCALT, AVEKAPHFXM, MV S/S/) 

EAMNMAT (IAVIT,IEVE, A) 

FAMMAT (IOUT FOANCY CANAM FACTOR, FAR RAIA) 

FOOMAT CT VeI Mel (AK oEI1A,4)) 

EQRYRT (CTH ePW bie gtehVolO(T Pel AVISIN AH o/s) 

FOQVAT (708 THE StIN AVED ALL Ve FOR Vemwg Ag AHI gFIAAAY 

FORMAT (79H THE SUM OVER ALL Vem FOR Vemol3e2HISoF 17 ef) 

FORYAT (2747 P CENTRPATA FACTORS (CANGSTPONMS) FOR gRAI0) 

FORMAT (Fuel rl IX oF hreh)) 

ENPAT (ott y 
EQRMAT (°st SISSY 
EANUAT (IFT PED CEMNTONING FAPZRAIAY 

ENQOMAT (90643 DOHA CENTONING EAP, RAIC) 


56 


P>?>P> PY PP Pr PPPF> PPP Pr > D> PPPPPPPPPPrPPrP PrP PrP PPP PrPrPPpPr vs, PrP > PRP > PP PP bPPrpP> PP 


LA 
497 
40808 
49a 
RAR 
aA 
and 
ana 
BAL 
aAK 
BAK 
aA? 
ann 
g&AO 
R})A 
®1) 
512 
813 
814 
&16 
S14 
817 
818 
819 
&2?0 
821 
&>2 
823 
R24 
R28 
&>6 
a>9 
GR 
@20 
ean 
Ra} 
q19 
R42 
R34 
R41e 
£46 
RA7 
gaa 
220 
&4N 
84) 
8&4? 
£44 
R46 
Bae 
£44 
247 
gue 
€49 
aqA 


118 
11° 
120 
121 
12? 
123 
124 
128 
126 
127 
128 
129 
139 
13) 
132 


133 


FORMAT 
FOPMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
TALEM,) 
FORMAT 
FN 


(79H] P#H4 CEFNTROIN FOO RAIN) 
(RAIN) 

(2344 FRANCK CAMAAN FACTORS } 

(aray 

CA07TAePAK oF IFS )) 

(13H © CFNTROIAS } 

CACIKeATAGFRAL)) 

(14H PHFD CEMTENINGE } 

(16H R#8a CENTENINS 

CLHON AX oe AHVVE eo WKole( 170 IX ISPKOLHV/) 
(24H1 PATIO OF CFNTPOINS FOR eAAIN) 
(THeavyelal*XeFhes) } 

CTH oP PetmKMol4lAVaFheh)) 

CTH eK TA( BK eK Keb)? 

(26H SCHR AAT SUCCESSFUL FOR I2e2AHTH TIMEs GO TO NEXT PPO 


(28HO SCHR DIN NOT CONVERGE FOR I2017HTH TIME WHEN IT = 13) 


157 


PPPPrrPrDrPrprYrppprpY»y»ppr»r p p 


65) 
55? 
Lam! 
TY) 
655 
Res, 
ea? 
a6a 
REO 
86N 
561 
56? 
B62 
R64 
&H5 
£64 
G47 
R6R- 


QIAROIITINE PATETT 


COMPAR KE(A DeYT(H4) XN ITIVIG) OVI IAINAN) OSI ITOANON) GETRIAL (AND OQTRA(2 
Lye DTEN CUD SITS CPD eZIMSE( A) eZITVA? } o 7 TEN( 3) oFCALC(30) oXPRN(5) eBVE (IN 
PreVPON (BK) OTERMTC IC) NHEMT CIS) oRVEC IN) eRESULT (39930) eRCMNTRD( 30549) 
FoRRCON( AO I3%) eRORCN( 300399) oRRRRN( 31030) KV (490) ONE 2) ONL (2) oDF (2) oPRE 
GMT CLO) PRET O19) DIV(50) oStIMPP (30) oSUIMP( 30) oTITLE(B) oXHoNSTAOMSTA oM 
BS eXMIN OXMAXSINSC oMBEG oMUISED ONT ONS oMAXIT oFAC% oO 2M 9 IKON SERS oI PSIQ0SCH 
BPeSUCTIINN) 


AIMENCTON ?» 
weArTayrttay 
WV myeeye 


(BOV—e AXLAND © PVIAMIG CKIAND» AKIKA) 


CHECK FNS FYTPAPOLATIA® AT SMALL Re 


TE UX TeX! 
WMAK SKIT (1) 

"SCe#) 

ILL ANC Tena) 
CALL XTeAOD 


Oat aX 


C199 Y920? 


(XP eV oNLL eK MoV amo XHe XMIN OXMAX oINSC oFACM oMRES ONE ( TKON) ) 


WV IMNEXKN (MT yp eved 


URE 2 **e} 


PREPAOE TO INTERPOLATE, 


ve 
YMA SXMANT 


TF (XKMAX=YTOMLE DD Goto 
WMAK SEY TONLE) 


CALL SPLYY2 


(XT eVT o7 NLL eo AY PY oC KONK) 


MEN Me MREGSINT( CXMAKRX“TM) /XH) 


YOOaKMye 


AH B& TT sMFG teen 


CALL SPNXY7 
vi ttlyeva 

WOCT PY 2¥Nn 
toler Basler sac) 


(XT oeVLe7eNLL o¥O00V0) 


CHECK FOP FYTPAPOLATIAN AT LAPRE Py 
TF (YMAVTOYMAN) 757K 


URE Gs Eanes) 


VMTN YN (MAF G=H) YOY 


YVAN aWMAXT 
INSC 2? 
CALL XTPAP 
WNC” 

wT Nearer T 
tag ttgunNrnrot 
DE Typ 

ENA 


CXT oe VT eMLL eX MeV oo YHo XIN OXMAX OINSC oFACM MARE” NE TKOM)) 
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DvDoasanoarsadvavwwdwInanatvpoaDPMPBIRABsadDvsswBABrPBIvVvoBsvPWMWPPBMIvwIBIDIwNPP 


ORDAP AE BNO 


1A 


SURROUTINE SPL YYZ (VeVe7eNafhehol Ny) 
ATMENETAMN Vo TPe VEDde 701) 
ATVENSTION AIT) 6 P1100 CEL) e® MYC) 
PATA DA/OITDEE DEDEDE DEDEDE DE TS 
MATA NA/NITFIKE PEPE ZEDEDR DR DE DS 
NP aha} 

nQaN=-? 

AO PD Kalen 

Nx(M pay (rel peren) 

TE (9Y(K ID Dele? 

PRINT Be YIM Ve VIM SH} PMV IMN Dg! 
RTAP ANY 

CONTINGE 

A(T) #-Nv(1) 

COLI s—Le/N¥( 2) 

RCL elLe/NK(TI=—C(1) 

707 eM, 

AY Me(CVIPAVHKYCLIIZAYID) 

AO 4 Kare MP 

vv ste} 

ACK) BNY (KK eNnK 

A(Kys (Ny (KsNW (KK) yen? 
C(KpeNy (Kr pane 

AYXOa hve 
NVKE(V(KO1PHRVOKDP/NX IK) 

(Kp ANVvehvyP 

CONTI MHE 

AAwHVe/AQKIMD) 

COeMKVe/AY(NMY 

MARe-( AASCC) 

ACAD RAAFA IMM) OAC NPY BAN 
AECNPBAASC(MDYSA CNPP ELC 

TUN) wAAHZ (MD) 

WE eA(PPHA(1) 

AC APRR(DPeKF HAT) 

CUP YBCO PAV—XFBC(Y) 

NN hk argh 

vUst4} 

YFeA(KK)/P (KW) 

A(YK YEAH KF EC(W) 

7UKK STUY JeVE HT) 

CONTINUE 

FOENPEZCM PSO CMY 

RQ & ve) ghn 

vet eNol 

7(KM ye(7 (HE eZ (KYO Er(KrK) Y/N (KH!) 
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APPENDIX IV 


ADDITIONAL TABLES FOR BIRGE-HOPFIELD TABLES 


Tue purpose of this appendix is to supplement data appearing in the main 
text of this technical report. New data has been published and these data have 
been incorporated into the tables of Franck-Condon factors included in this 
report. To minimize the changes to the main body of the report, these data 
are being included in this appendix. 


The confusing and contradictory tangle of excited states for the very high 
energies in the nitrogen molecule has been unraveled (Refs. 41, 42, 43). Dress- 
ler (Ref. 41) reports that the former twenty states of excited nitrogen can be 
regrouped into three valence states and three Rydberg states. The three valence 
states are the bins bly, and d'. It will be noted that the first two states 
are those associated with the Birge-Hopfield band systems which are pertinent 
to this report. The renumbering of the levels of these states has vastly 
increased the number of vibrational levels for the upper states of each of these 
systems. This necessitated recalculating the vibrational transition probabilities. 
The input data used to calculate the vibrational potential for the b'!r state 
was the deperturbed spectroscopic constants of Dressler. These data are listed 
in table 121. For the input data needed for the b/n state, the data given by 
Carroll and Collins was graphically smoothed and used. These smoothed data are’ 
listed in table 122. The potentials resulting from these data are listed in 
tables 123 and 124, The potentials were generated using program TURNGPT. The 
vibrational transition probabilities were calculated using program FRANKON. 

The Franck-Condon factors are listed in tables 60 and 61. 
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Table 121 


SPECTROSCOPIC CONSTANTS FOR THE Nob" rt STATE 


vy Sy *y 
0 371.9 1.151 . 
1 1112.0 1.143 
2 1848.0 1.135 
3 2580.0 1.127 
4 3308.0 1.119 
5 4031.0 1.111 
6 4750.0 1.103 
7 5464.0 1.095 
8 6174.0 1.087 
9 6879.0 1.079 
10 7580.0 1.071 
11 8276.0 1.063 
12 8960.0 1.055 
13 9632.0 1.047 
14 10292.0 1.039 
15 10940.0 1.031 
16 11576.0 1.023 
a7 12200.0 1.015 
18 12812.0 1.007 
19 13412.0 0.999 
20 14099.0 0.991 
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z 


POTENTIAL ENERGY FOR THE N>b' id STATE 


V(cm~!) 
371.9 
1112.0 
1848.0 
2580.0 
3308.0 
4031.0 
4750.0 
5464.0 
6174.0 
6879.0 
7580.0 
8276.0 
8960.0 
9632.0 
10292.0 
10940.0 
11576.0 
12200.0 
12812.0 
13412.0 
14099 .0 


Table 122 


Tain 


1.3672 
1.3158 
1.2825 
1.2568 
1.2355 
1.2170 
1,2007 
1.1861 
1.1729 
1.1607 
1.1495 
1.1389 
1.1283 
1.1181 
1.1083 
1.0988 
1.0896 
1.0807 
1.0719 
1.0640 
1.0701 
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zt 
max 


1.5282 
1.5953 
1.6443 
1.6859 
1.7231 
1.7576 
1.7899 
1.8207 
1.8501 
1.8786 
1.9061 
1.9333 
1.9612 
1.9889 
2.0167 
2.0445 
2.0723 
2.1001 
2.1283 
2.1552 


2.1504 


Vev) 
0.0461 
0.1380 


o 0.2292 


0.3199 
0.4102 
0.4999 
0.5890 
0.6775 
0.7655 
0.8529 
0.9398 
1.0261 
1.1109 
1.1942 
1.2761 
1.3564 
1.4352 
1.5126 
1.5885 
1.6629 
1.7480 
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Table 123 


POTENTIAL ENERGY FOR THE Nob! n STATE 


y Vem!) Fin Fnax = (ev) 

0 313.2 1.2053 1.3819 0.0388 

1 947.0 1.1553 1.4540 0.1174 
2 1632.0 1.1288 1.5031 - 0.2024 

3 2354.0 1.1090 1.5445 0.2918 

4 3110.0 1.0969 1.5792 0.3856 

5 3889.0 1.0859 1.6121 0.4822 

6 4682.0 1.0765 1.6437 0.5805 

7 5481.0 1.0687 1.6751 0.6794 

8 6278.0 1.0616 1.7063 0.7784 

9 7063.0 1.0550 1.7390 0.8757 
10 7831.0 1.0490 1.7725 0.9708 
ll 8578.0 1.0432 1.8070 1.0635 

12 9303.0 1.0382 1.8426 1.1533 , 
13 10004.0 1.0337 1.8794 1.2403 

14 10682 .0 1.0298 1.9169 1.3243 

15 11337.0 1.0263 1.9555 1.4055 

16 11967.0 1.0231 1.9954 1.4836 

7 12575.0 1.0206 2.0359 1.5590 

18 13159 .0 1.0182 2.0780 1.6314 

19 13719.0 1.0162 2.1216 1.7009 
20 14255.0 1.0150 2.1672 1.7673 
21 14768.0 1.0147 2.2149 1.8309 , 
22 15257.0 1.0135 2.2652 1.8915 

208 


AFWL-TR- 70-108 


Table 124 


SPECTROSCOPIC CONSTANTS FOR THE Ngb?n STATE 


Ie 


an NT BS 


Sv 


313.2 
947.0 
1632.0 
2354.0 
3110.0 
3889.0 
4682.0 
5481.0 
6278.0 
7063.0 
7831.0 
8578.0 
9303.0 
10004.0 
10682.0 
11337.0 
11967.0 
12575.0 
13159.0 
13719.0 
14255.0 
14768.0 
15257.0 
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